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'TOL IN GROUND EFFECT FLOWS FOR 
CLOSELif SPACED JETS 


David Migdal, William G. Hill, Jr. . Richard C. Jenkins 
and Michael J. Siclaii 

Grumman Aerospace Coi’poratlon 
Bethpage, New York 

1. SUMMARY 

The primary purpose of this study was to obtain detailed pressure and veloc- 
ity field data for twin jet configurations In-gi'ound -effect and to develop flow mofiels 
to aid in predicting pressures and iipwash forces on aircraft surfaces. For the basic 
experiments 50.3 mm (2 ini diameter jets were useil, oriented normal to a simulated 
ground plane, with pressurized air providing a jet velocity up to 90 m/sec. The ex- 
perimental data consisted of (1) the effect of jet spacing and height, on the ground and 
u'jwash pressures (2) oil flow visualization photographs to identify streamline directions 
(3) the effect of simulated aircraft surfaces on the isolated flow field (4) the effects ot 
adding a third ict on tlie two-jet upwash properties arzi (5) the )et irduced forces on a 
scries of two and three dimensional botlies with strokes and variations in fuselage cross 
section. 

The isolated twin jot ti'sts revealed regions where the lonntain flow was un- 
stable, I'eUnv ambient piossuie rogdons n the upw.ish, and the formation of giound 
stagn.ition lines (.ind upw.isli How ) even when tlie outer boicidarios of tlie free jets have 
merged. With the simulated tuselage surfaces jilaced Ix'tvveen the icts significant lift 
fences weie obtained in-ground-cffcct with the addition of strakf.s and squaring the 
lower fusel. igo contour. Similar eftVets were not otiLuned with flat plate moclels. 

A now flow motiel was dcveloix-d for \ eriL-allv oiiented twin jets whicli ac- 
counts lor the non-fullv foi mod wall lets .it clo^f sjxicing. It includes a jet imiiingc- 
ment /one w.d! )ct ti.msition uuxicl, tiie iffcvis of fuselage cross section on fountain 
induced UU I'oiccs, and .idequ.dch picduts the elteets of sjuicmg and heiglit above 
mound on Row nelil rressures and npw.isli toie»".. 

It IS ceiicluled tl'.nt (V) flow iield m.-taluluies enn present at low heights 
a auc giou'id: (2) in a Ivin ;ei com .giiranon the onset of |Cl boinnlarv -ncrgit'g docs not 


preclude the formation of fount.*ins; (3) values over 207c lift/jet thrust can be obtained 
in-proiind-effect due to foimtains impinging on the underside of a simulated fuselage 
for a twin jet configuration, depending on the shape of the fuselage and dimensions of 
horizontal strakes and; (4) the use of a modified wall jet transition model adequately 
predicts the trends with height above ground and jet spacing on the ground and upwash 
pressures. 
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2. INTRODUCTION 


Lift and control for V/STOI< aircraft operating In-ground-effect presents a 
cidtlcal condition in sizing the pi-opulslon system. Ground pro;clmlty effects can re- 
sult In ijvre jet Induced lift losses, or produce positive fountain lift, depending on the 
aerepre,.'!' ion configuration. The complexity of the I'csulting flow field and the sen- 
F’tr/ity to many design parameters gives rise to a large body of exiierirnental data 
(Reference 1) which model various aspects of the flow field. References 2 to 6 deccribe 
the development of one technique to predict the pressure distributions and aciodynamlc 
characteristics based on a modular approach to the various flow regions shown in Fig- 
ure 2-1. The four major subregions consist of a free jet, impingement zone, wall jet 
and upwash flow. As shown in Figure 2-2, all of these major elements arc present in 
the two jet ground Impingement pi-oblcm and the basic two- jet sc ution can thus be 
used as a building block for multl-jet configurations. 

Tlic simplified modular approach utilizes Ixisic mass, momentum and ener^O’ 
con'^enation principles along with empix'icai data desenbing entrainment of ambient 
air (c.g. maximum velocity decay and velextity profiles as described in References 
2 to }), A review of the baste experimental data regarding two-jet behavior iixiicates 
several limitations when the let spacing is decreased. .■\s reported in Reference 8, 
instabilities in tlio upwash ilow weie noted for some combinations of spacing and 
height. .Also, as the jets arc’ brouglit closer an upwash is formed prior to tiie csti.b- 
Ushmvnt of a fullv de\olojiod wall jet. For even clof'Cr spacing^, tne upwasli can 
vompletelv disappear when the two free jots merge a.s tliev grow in size due to en- 
trainment of ambient air. The object of tins study was to cxjiloic in more detail, tlic 
tw’o-ict close spacing prolilems experimentally to help establish a more accurate model 
of the resulting flow field. 
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3. SY^lBOLS 


b 

AB 

d , 1; = 

6 (deb^) = 


AF 
h, H 
L 

M, M 

Pt 

Ap 

P 


Boundary radius 

Body depth from nozzle exit 

Nozzle diameter 
Boundary layer lieiRht 
Interference Force 

Jet height above ground 

Body length 
momentum 
Total pressure 
p - pa 

Static pressure 


Ap 

Q, q 

p (rho) 


P - pa 
P . - pa 

O 

D>Tianuc pressure 

Density 

Radius 


S = Spacing beb-vc-en jet centerlines 

a -• Strakc depth 


T 


Jet tin ust 


\' 

W 

W' 

c. 




\'eiucu\ 

Jit exit %elocity 
i’ody uidth 

Distance betv ecn strakes 

Distance perpendicular to line j.xning 

Distance along line joining jet centers 


- P.elic‘ height aho\ e ground 
= Dista.ef ,• frijpi ict exit 


;j-i 


jet centers 




Z 

Z’ 



SUBSCRIPTS 


a = ambient 


cs - stagnation line center conditions 

g = ground plane 


j.J 

JH 

m 

i\ 

PC 

s 

0 (zero) 

.5 

.25 


= jot exit corelltions 
= jet half-width 
= ma>:imum 
= nozzle 
= potent! al core 

= static 

= initial wall jet radius 

= dimension at 1'2 maximum pressuic 

= dimension at l/4 maximum pressure 
(1/2 maximum \eIoclty) 










4. MODELS AND APPARATUS 


The experimental portions of this stidy were performed In the Grumman 
Research Department- Two similar air flow facilities were used, one for flow surveys 
(desigrialcsd Facility 2) the other for vehiele force and pressure measurements 
(Facllitn 1). Sketches and piiolographs of Uieso two facilities are sho\>.'n In Figures 4-1 
througl -9. These facilities utilize centrifugal fans which are belt drhen by 220 V/3 
phase "--1/2 'IP motors. The flow rate and hence total pressure in the settling chamber 
is cont< oiled by throttling tiic inlet to the fan. Experiments are usually conducted with 
a tola'' pressure minus atmospheric pressure of 5004 « a fiO. 1 inches of water) result- 
ing in an exit velocity of 91 m/sec (300 ft/scc). 

The flow from the fan into the settling chamber is non-uniform with iiigher 
velocities lowaids the outside of the fan (top). This flow Is passed through a short 
constant area section where vanes redirect the flow for a moie uniform distribution. 

The fan is connected to this section with a flexible material to avoid transmission of 
vibration from the fan to the settling chamber. The air then flows through a diffuser 
into the Large settling chamber. To br»- j the length of the facility within the available 
space, large-angle segmented diffusers are used. These diffusers are an extension 
to the tw’O mmcnsional work reported by Kline (Ref. 7). The diffu-ser of Facility 1 is 
upproximatclj n 33° square pyramid filled for 2/3 of Us length with 7“ paper cones. 

The flow at Us exit is almost ooiiipletely attached and uniform. 

The diffiLser of Facility 2 b'. a 00° pjTamid and filled with pjTamiclal sections 
of plexiglass (sec Fig. 4-'5) of approximately 7° angle. Flow into the settling chamber 
is liion coiitlitioncd by honeycomb and screens and passes into the settlbig chuimbcr. 

For tlio pairs of 30.8 mm (2 inch) diameter jets in most ol the current experiments, a 
t>10 mni (2 footl square cross section leads to an area ratio of 32 (Facility 1) while 
raeility 2 has an area ratio ot 2''' This produces a very uniform flow of low turbu- 
lotice le\c! (less than l/2''o Rf.iS). 


To facilitate cnan.^iiig the nozzle spacing, a p^iir of nozzles is mounted ui.c,ori- 
tiieall'. on a pair oi liu-cs. whicii are clampod to the fionl face (r’gare 4-(.), providing 
eonlinuous \arnlion of nozzle spaeing. The nozz.les ,>re spun aluminum wifii an cn- 
tiante contour identicalto an AbME long radius fie'' metering iio.''zlc, but follow eo bv 
1 cn'i-^tant area sectiO’'> 3.3 diameters lorg. .'ihgument of nozzles in horizonai duce- 
:.on IS done v.it.a icfcioimr to gra'.ity bydrop line and buublc level, and wnh reference 
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to oilier Jipparntus (traverse etiuipment prln’ailly) by cross hair sl>;htlngs. A chcv^k 
on tlic allgiimcnt Is porforn'iecl by tra\ersln)i tlie free jets and tracklnf; centerline 
(maximum ci> locations usln^: an electrically powered variable speed traverslni^ mech- 
anism. 

Model s 

The m • lc!s used for the upwash body Impingement expei'lmcnts Include sev- 
eral clasi'os. The Initial experiments on the effect of jet spacing were conducted 
iislfig the two-je; configuration shown in Figure 4-7. Cylindrical fuselage models 
(Figure 4-S) with a range of corner radii were iiscxl for the effect of body shape. 

These liotlics wire also Instiumcntt'd for pi'ossurc distributions and were tested for 
\nriou.s htixike heights. In addition to the basic tests. The strakes are .«hown on a rec- 
tangular body (I-'igure -l-D) whoi’e the nozzles weix' eontilned witliin tlie planform, A 
third elans of IkkIj Is a large tl21) x 12D) flat plate. Dimensional daui for the \arious 
lodies aie eoninmod m Figure 4-10. 

In - ti unu.T.tation 

'i lie coouiin.ae system for locating no/,/lcs md pixilx s Is sliowti in i’lgux’e 
1-11. '1 he position of the giowui nl.me and of Uio ptooes ai-c mea navel by the voltage' 
oalput of a 10-tnrn potenlio'iu'ter loi most dati in tills ’-ciioit itmicss othensiso notCvi). 
For !iieas';re'iiu'iit of siiu'.U di meiusions’, such as the I jtind.ay la\cr thickness, 
r'l' isurements were made with i linear pob-ntiometer cr a vernier cali]X'r, For the 
ti ueise nnils, with screw tliread elrive or rrek aixi pinion drive, (all probe movements 
.md the gtoiuvJ plane motion on Faclllt> 1) a rotary pot l.s used to measure cable motion 
t oii!'f..ted to ilie ino\-aig .ippai-atus. For the gioond plane of Facility 2, a rot.UT pot 
was .itt.ielied to .i .sprocket .it one end of the chain loop, hiiicc some hyslcreM;- occurs 
for iiu of liu’sc systems, ti'ivorses were taken onlj' in one direction, and a voltmeter 
0 .'V u-^ed (.hc-ek .m Iml'.il output at .i starilng location and jicnotlicalh to cheUv 
. .!i>'i .ilion. 

i h.e puipio'-i' of pio', iding an eb'itric.il out^nit for position is to proouec a ton- 
tiruous I'un.e of tlie p.ii.nmti'r l>eiiu: nie.tsiui'd ('ehkle forces, jMcssurL'', )’ot wire 
\i ioiitu'.s, lie,) Noisus po'-liion, .mti lo plot ihal euiwe on an X-1' Uccoixler, leal- 
liii.L aiuilog d itii icduttion tccl'incpie was used to pixstiut ag.imsl missing uiuuuial be- 
h.uioi (,if tre (l.il.i oi the dc', eloi'i-icnt ot piol'K-ms diiiing a t< st. 
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rresf.ures were iiieasurod with differential pi*essurt* transducer. For the 

pltot-statlc probe, the total and static pressui*e tubes wei'e each connected to separate 

tr-nsducor referenced to the atmosphere (an ureilsturbed area of the room). Fach of 

these pressure tubes also contained a "T", and a third transducer was connected to both 

pressure.'s, thereby directly reading P. - P or q (another analog data reduction toch- 

^ 8 

nlquo). The pressure transducers re calibrated and periodically checked by use of an 
oil manometer: vertical or inclined, depending on the transducer range. 

Vehicle force moaauremente were made using a flv^e -component strain gage 
balance, while force measu.'emcnt on the groiuid plane utilized three separate force 
elements, spaceKi at a distance from the force center. Calibration of the force bal- 
ance was achieved by placing weights from a laboratory scale, set on a special gradu- 
ated b(Kiy, and spot checks w'ere done periodically on the model being tested. 

During efforts to Identify the cause of the measured laige nejcatlve static 
pressure differentials, \olocities ueix.* measured witli a hot film system. This .-lys- 
tem Is essentially the same as a hot >%tre system; out measurem ents Involved a con- 
•'tant tempoi'uture .inemometer, a slgn.\l llnearlzer, and a true RMS tune-avorripng 
\oUni('ter, Output from t!ic traverses a as .dotted on the X-Y rccoitier. 

Oil flows woix.’ used for loth ground plane Indications of the end of upwash 
tonnatlon \vlt!i Incnnising jet interfeixmce and for qualltatlxe flo^v description \^lUl a 
splilu r plate. The oil Ivlso uas 10U'-;iO motor oil and it v. as coloied for Idemiflea- 
tlon and tlilckened by .utist oil colors. 
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5. JET PROPERTIES 


Figures 5-1 and 5-2 contain the results of traversing the jet plumes with pitot 
pressure probes at various distances (Z') downstream of the nozzle exit plane, in the 
absence of a ground plane. The cliaracteristics of the initial jet can be inferred from 
tlie data very close to the nozzle exit (Z' /D = 1/16). As shown in Figures 5-1 and 5-2, 
a relatively uniform (top hat) profile is available at the jet exit with the eccentric jet 
spacing mechanism and tapered settling chamber. The Figure 5-1 data Indicates that 
for S/D = 2 the jets start to merge 4 diameters from the nozzle exit and that the 
merging is nearly complete at 20 diameters. Merging is considered to be complete 
when there are no distinct peaks in the profiles. These peaks distinguish the central 
region for each free jet. Thus while the jet outer boundaries spread at an equivalent 
angle B (tan B = (S/D - 1) /(2Z’/D)], and for S/D = 2, Z’/D =4, B is approximately 
7®, the inner core region for each jet persists for a considei'able distance down- 
stream. 

Static Pressure Inside Nozzle Near Exit 

Static Pressures were measured inside of a nozzle at Z'/D = -1/16 (up- 
stream of the exit) to investigate the effect of the upwash on nozzle exit conditions. 

Two px'cssure taps were used, located on opposite sides of the nozzle. The static 
pressure dilferential was less than 2% qj for all cases. 

Figure 5-3 (a) sliow's the variation of nozzle exit pressui'e as the ground plane 
height was changed. The smooth curve shows the variation with grotuid height for 
single jet impingement, illustrating the influence of ground blockage when h/d< 1. 5. 

The other two cuires for dual jet impingement show that the effect of the presence of 
the gi'ound plane on the nozzle exit nmA' are noticeable for h/d< 4.5 Data taken on 
the inner (towards the upw’ash) side of the noz.zle illustrates an increasing pressure 
as the gi'ound plane gets closer, which appears to be caused by interference betw’een 
the incident flow and the upw'ash, Tlie pressure rise with decreasing h/d is much less 
for slighily larger spacing (Figure 5-3(b)) and disappears for vvide spacing (Figure 5-3(( 

Foi a given nozzle spacing, decreasing h,/d will decrease interference iue- 
iwcen upwash and incident flow, and at some value of h/d such interference should 


disappear. Static pressure surveys showed that, in the absence of interference be- 
t^veen upwash and incident flow, the static pr 'ssure on the ujw.ish (inner) side of tlie 
nozzle should be slightly more negative than on the opposite side. Figure 5-3 ^a) 
shows that the inner tap pressure drops I'elow outer tap pressure when h/cl = 1. 75, 
which should represent the ground height below which no Interfei-ence between upwash 
and incident flow exists. Figures 5-3(b) and (c) show this \ due of h/d to be 3 for 
S/d = 3, and 4 for S/d - 5, 6, 



















G. GROUND TLANK PnOPP:RTIE:S 

ITie upwash formation rcRJon on tlie ground was Investigated with oli flow 
pnttczvs to aid In the identification of conditions uOiere t)ie upwash would no longer 
form l>ec;ause of jet/upwasli Interfci'cnce. Several concentric rings of oil dots were 
placed around each jet Impingement j>oint and the flow was allowed to run for several 
minutes to develop a stable pattern. Figure 6-l(a) shows a tjidoal radial flow from 
the jet Implngenient points meeting to ^orm a straight stagnation line. For many 
cases, a double line was found. F urther study showed this to lx; caused by a separ- 
ation point. Oil drops placed wltliin tins bubble ran outwaixls along tlie outer line and 
towai'ds the two oil lines (Figure G-l(ld)- Insufficient data exists at this time to dc- 
ternune the flow conditions ttut govei'n tl e formation of this separation bubble. 

Figure G-1 shows a t erics of oil Rows for a spacing of S/D = 2. vO (closest 
araliable from ccccntiic no7/5les). Starling at a height above {jround of li/D •- 1.5, a 
patter-n of two separ-ate radial spreading regions with a central stagnation region oc- 
cur’s until an H/D of a!.out 20. This is generally a!x>vc the region of interest for rd- 
feeting aiqdanc Interference effects. Note tliat the separation region disappears at 
about ll/d ■- ’.0 (Figure G-ldi). 

Stream sur\e,\s w itiiout tlie ground plane, taken by traversing a pitot prol.K 
across both jet centcilincs, snowed tiiat the two jets had pailinlly merged 10 diame- 
ters from the exit. Grcuaid plane prcssui'c profiles at II/D = 10 mdlcatcd a pressure 
rise coi res pond uig to the stagnation line. 

For li/D = 15, \%!ieiv the free stream suivey showed two )>oaks, the ground 
oil flow pattern show oil two <ilstlmt jet impingement points and a sUignation line also. , 
At li/D = 20, where tlie five stivam pitot suiwej showed no distinct jK-aks, the groiuid 
oil patrem also Indicated tlmt the two jets had nicrgtd. 

It t’lerefore seems leas'on.ible to use the rule of iluimhthat: if separate 
pt a!;s exist in t!’c jntot presume profiles of tee free lets wiihojt a ground iihino, sep- 
arate impingement and sta'inauoi! line/eowash furnustion will occur on tne ground 
plane. 
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Groiunl Plruic I’rossiiro'-; 

Ground prossuivs were obtained witli a line of pressure taps pcipendicular 
to the line connecting t!ic jets and sj^aced one-half a jet diameter apart (alonjC the k 
axis of Flfjurcs -l-ll and (i-2). The t;aound plane was mcned slowly in a direction 
pcipendicular to the line of taps and to the nozzle axis producln}? a continuous display 
of groimd pressures. Tiie ^^round cooixllnates x and y are oriented parallel to and 
peipcndloular to the staj^natlon line, with the oiigia at tl>e upwasl>. center, Fijturc 
6-3 contains a complete set of fp-ound pressure data o'ltalncd wltli various jet spaelnp:s 
and hel(;hts abo\’e proiu’d. Fipivc G-3(a) contains the ground pressui'cs for botli the 
slii(;le and dual jets at S/D - 2, I'/l") “ 3 and demonstrates that an upw'ash Is formed 
(peak pressui'e at Y/D ~ 0) prtov to the end of the lmpln):ement zone. The data im- 
plies an ir.lioaixl sWft of pressure p/eak and a distoiilon from axial sj mmetry when 
the second jot is turned on, 

I'itxurc G-.!(li) contains a complete data set for the closest spacing (S/A") = 2) 
and ai vlnc, noz/le l>eit;ht abo\e ground. 'I'he general trend Is' that tlie peak \ahios 
under tae jot and at the stagnation point (V/D - 0) decicasos wltli Increasing h/D 
while Dk piotile i'roadons. Nolo ihat laiiler the jet (\ /D ■ 1, Figure (;-3(l))) t..e im- 
pingement ino'-suie Is \\iiiun of the jet stapiation picssure for \alues of li/lJ up to 
1, As h,^) iiu'i'cases tJie jet impi'igemcnt pressure decn\s in a manner similar to the 
fiei' u'l <iec'iv in ma\imum pres'^uie (i'limre 5-1). In gencrul, the maximum wail 
pitssuie IS tiigher Uian the corresponding \a.lue for tlie free, jet at the same distance 
fiom the nu//l(' exit, I'lgiue !’>-3(e) contains me ground plane pressures for s/D = 2.5 
;mk 1 a laiu’a' ot iio/zle li/D troni 1. 5 to 15. Li tins, as well as otiier eases to follow, 
sultleient da, a i> oroMCied to determine (lie mmximuin jn'essiiro on tlie center of the 
stagiiation line ^ (r'giiro 6-2). Note iliai increases with, ciecrcasing h/i) until 
a!' !i I) of 3.0 the -.Inm in 1 igiire 6-3(a) us similar in tliat the laiue of P dot's not 
eont.mi ,)1\ nu le i-m.- ■ ith dec I'ca'-im; li/lk In addition at 1;/T) - 2 (I'lgiue 6-3((l)) it 
w.i" net ros-'iule to ol tda a -.taole profile .larl two pt'.<k '.alae.-. weu' cxnlbitofi. 'i'liis 
ii- it il'illi’, V, no: e'v oix ed at other ■ p-wings. At .''/H - 3 and ii/i> = 1 . .5 to 2. 0 tiw 

u;a .ml, wo^'h! mgv- fi.-.n. cie’ :^ilU to tSf ctlicr with no .ipp. treat distui ij.mee up to 
..'iH'ut L’a '-c I'l'iid'^, It eoalc' ix' so. itel'ed iici’"Oiie side to tiie c-ther with a tempor.irv 
b'uukage lU'Ui ihe iii'i/,'le c \u pi. me. 



This unstable behavloui- at SA> = 3. 0 Is shown In Figure G-3(d) with the ac- 
tual pressure trace at h/D = 2. 0. The large oscillations help to explain why in this 
case thei'c is a drop In the maximum stagnation line pressure. Since the probe Is 
reading a time-averaged value, and for a portion of the time Uiere is not jet flow over 
the probe, the time-averageu value is lower. 

Figures G-3(e), (f) and (g) show the grouiid pressure variation with Increased 
spacing. Note that as the spacing Is Increased (e. g. , Figure 6-3(0 for SAl = 4. 0} 
the]‘e Is less effect of the fountain upwash (at Y/D = 0) on the Implr^gemcnt zone 
pressure distribution. 

Summary Data of Ma ximum Ground Pressures 

Figure G-2 illustrates the ground pressure profile along a line passing 
througii tiie two jet Impact points and the center line of the upwash. The maximum 
pressure in the jct impingement r-egion was always greater than the maximum 
pressure on the center ot t!ie stagnation line (P^,, ). The dimension y . I'epresents 
tliu lialf w Idth of the pix:ssure profile at half the maximum stagnation line pressure. 
Along the stagnation line the maximum ground pressure (P ) decreases with distance 
(\) from the ccntoiiinc. 

Figure G-4 shows tlie filiation of centerline stagnation pressure wLtli nozzle 
spacing for a nozzle height of three diameters above gi'ound. The cur-ve shcAn for 
companson was predicted by assiuning that the ground pressure on the center of the 
stagnation line is equal ' * the maximum pivssurc that would ex'ist In a wall jet piro- 
file at the same rxidial location from one of the jets If the opposite jet were turned off. 
Tlie predicted ciiiwe is given by tlie relation 

p _ p 

pc, T .1 

— i. = 13, 7S(S/D) ‘ (1) 

- P 

^3 a 

V. lu-re P is Die stagnation pi'essure at the nozzle exit and P is ambient pressure. 

J a ' 

The eonstaiTs in tins cciuation weie deiired from pitot probe measurements taken in 
the wail jet formed b\ impingement of a single 4 inch diamcvcr jet. I'lguro 0—1 shows 
that this anal^si'. prosides an aclccjuatc jn'ediction of the contorhno ground pressure 
for jet spaeings of about 4 dianictcis or grearer. .-\t closer no/.zie sp icings the nieus- 
uiod Millies fall Irclon the picdietid cutae. This dev.alioa should lx; exirceted because 


i>-3 



at S/D = -1 the staRTiation line Is two diameters from each jet Impingement point, aiid 
this radial distance from a slnRle Impinging jet repi-csents tiic hoglnnlng of the felly 
formed wall jet region. Hence tlic measured and predicted values oegin to diverge, 
IxjcaiLSC at close jet spaclngs the gi'omid flow approaching the stagnation line Is not a 
fully foimed wall jet. A modified wall jet anal}-sls Is needed to pro\1de a better match. 

For a given nozzle spacing, the gi'ound pressure on the center of the stagna- 
tion line (P - P ) vailed with nozzle height altove ground as shown In Figure G-5. 
cs it 

For a nozzle spacing of G diameters above ground. Increasing nozzle height above 
ground brought about a gnsdual decrease In -(P ^ , - P ) tliat appears to be caused by 
Interference between the edges of tl>e ujwsash and the shear layers surrounding the In- 
cident jets close to the nozzle oxll plane. At closer nozzle spaclngs, this \nrlatiou of 
(P^^^ - P^^) \. itli height alKJve ground was more exiixune but followed (ho same trend. 
For nozzle spaclngs Udow G diameters, tlie offtets < . incident jet-upv asli Intel ference 
se\erol> dlstoited the entire gix)uiul plane pix'ssure profile at large ll/D (sec Figure 
G-G(d)). 

For each nozzle spacing die mK-.umim value of (F - ) occurrctl willi II/D 

CS 

between 2 aixl 4. I'lie decrensc in - P^) at lower nozzle height appears to indicate 
the onset of mstabihtv condiiionj and niai be caused by the presence of die grountl af- 
fecting nozzle CMi coiibtions. A severe instabilitv was cneounicred at a nozzle spacing 
of 3 diameters for ll/D loss tlian 'J-l/J, The grouixl piessurc profile appeaietl to be 
temporarily staldc in either of two unsymmetne shapes. This instability appeared tO' 
be peculiar to ict spacing of 3 dinmoters. For i irger or smaller ict spacing, the 
ground plane pressure profile did not o<hi!m this bistable shape for nozzle heights down 
to I- 1/2 diameters . 

Tile v.in.ilion of nuotinunn grourid plaiic piossuio - P_^) along tlic stagna- 
tion liiie Is siiown in Figure G-G, Data t.iken at ll/D -■ 3 at each jet spacing have tx'cn 
not mahzed i>\ tlu' corresponding eentciline vtdue and plotted vcrsits X/D to illustrate 
(ae V. illation of (Ills profile willi jet spaitng. When plotted vcistus X/S (Figtae 6-7) 

(lie ciatti toi oaeli nuz.'le siiaeinu tall close to a jiioulc wlilcli ctui tie fiUod b' ; 



. 1 - 1 


Eq. (2) was derived by assuming that on the stagnation line only the components of wall 
jet velocity normal to the stagnation line contribute to the stagnation line pressure. 

The constants In Eq. (2) were obtained from wall jet measurements with a single im- 
pinging jet. 

Figure 6-8 shows the sisape of the ground pressure profiles taken across the 
center of the stagnation line for three nozzle heights above ground with a nozzle spac- 
ing of 4 diameters. For each ground height the pressures were normalized by the 
maximum value on the centerline and the distance y was normalized by the half- 
pressure dimension (y -V The data show that the profile shapes are essentially the 
same and can be represented by a Gaussian profile. 

Integration of the ground pressure profile across the stagnation line provides 
a measure of the force exerted on the ground by the upwash, which can be related to 
the momentum £lu>; of the iqwvash at gi'ound level, Assmne that the gi'ound flow ap- 
proaching the stagnation luio from each of the impinging jets is an axisymmetric, ra- 
dially expanding flow that exliibits constant momentum per radian in the wall jet. If 
flow momentimi Is conserv'cd in the upwash formation region, the momentum flux per 
radian of the upwash flow normal to the ground at the stagnation line can be found by 
integi-ating the stagnation line pressure profile. Using a Gaussian equation to repre- 
sent t'us profile, Integration yields 
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Figure G-9 shows the variation of maximum stagnation line pi’essure and half 
pressure width with nozzle height for .SA'J = 4. The maximum ground pressure de- 
creases and tte width increases as the ground height increases. Tills general trend 
appeared for all jet spacings. Howeicr, as seen in Figure 6-10, the centerline up- 
wash niomentuni flax pei ladian computed from Eq, (3) using the data in Figore 6-9 is 
almost independent of nozzle iieight aboie ground, lliis result is impoitant because it 
coiiTums one of the simpiif> ing assumptions used in the prediction techniques. 


'I he dashed iint in Figure 6-10 represents tnicc the ilinist per radian of one 
of tnc incident jets. Jet Lhiust wa.s computed from velocity profiles that w'cre taken at 


c-r, 


the nozzle exit plane. We conclude tliat flow momentum Is conseiwed in tlie upuash 
formation region, and that the uiv.vash momentum at ground level cati be computed 
from the nozzle exit conditions. 
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7. UPWASII PROT’CHTIES 


Oil flow tcchnlqui’S wcie itsod to qualitatively assc*-'*! upwash flow behavior. 

A series of "splitter plates", (thin sheets of metal at right angles to tlic gx-ound plane 
that split, or divide the flow) %\cre used. In one sei'ics, the splitter plates were 
placed Ix'tueen the tuo jets along tlie centeiline of the upwash (In the X-Z plane) 
looking at the radial spreading chax-acterlstics of the floi'^. in t!.c second scries, the 
splitter plate was placed along the line connecting the two noz7,le centerlines (Y-Z 
plane) splitting each nozzle flow In hall, 

(11 1 Flow Results 

Tlic presence of tlie splidcc nlate will affect tiie upwash changing tlie flow 
Lioth liy damping tiie turhiilont fluctuations normal to the flow dli-cclion and Isy devel- 
oping a lioundary ia.er of lower energy flow which can Intcrfcxv when pressui'c gx'aril- 
ents aie cncouiUert''i. Wo tlieiefoic sought piimarily to look for the xx'gions of jet/ 
upwash xnterferonce with cietailod probe sun'oys to follow. 

For the splitter plate placed ivt\%een the jets (X-Z plane), and a wide spacing 
(S/l) = (’•), a uniform spicading flow was found at heights of H/O -- 2 and G (rigiire 
7-1 (li)). For nukh (.losi.-r spacings, ticglmung wicli S^D - 2 ue find a similar situation 
at low height'-, .\t Il/D - 2, (Figure 7-i(c)), the flow spreading ngx-ces well witlx oar 
radial flow moiiel, A1 a heiglit of li/I) - G, streamiines along Die center appear to lx; 
drawn in as if Iw entrainment fioin tin* jet (Figure 7-l(d)). At H/D = 3, a definite di- 
■ I'lgenee of tiie center stieamiine occurs (7-l(a)) anci going a little higlicr (H/D - ?, 

' Igurc 7-1(0) a definiU’ stagnation region is found whci'c tlie oil in the upper center rc- 
I ion IS flowing ciow'nwnrd. With another increase in height (tl/l) - 10. Figure 7-l(g)), ti 
'ociUion of tins stagnation point appears to remain the same, but the extent of Us effect 
i‘,c upwash IS liroadenod 

With tlie splitter jil.ite/iio7/ie a: i.ingement, we were aiili' to go to a spacing 
>" •'''D 1, i.e. , m''.7le edgi'v toiiciiing. At thu-> sjiaung, si aular effect is found at 

luwi, 1 i.ciglits, -\1 ii 'Vi - (1 7-l(h0, a -'n'vadlii'.. in the uppe i -i entral I'cgion is 

noun and a i.le u -tugnaiiun ie';i*.r .-ei-n .it l!/!> - 1 (I'lguie 7-l(i)), 




For the second set of splitter plate experiments (Plato in Y-Z plane), we 
again began at a wide spacing of S/D = G, Results are shown in Figures 7-1 (j), (k), 

(1), t>eginnlng at an U/D of 2 and increasing to G. At the height of H/D = G, the d>Tiapi~ 
Ic pressuie was so low that it was difficult to form the oil patterns. Note the large 
angle (about 30“) that appears to lx; tlic outer wall jet growlh angle. We believe that 
this is due to the comer flow' and is a warning on the qualitative accui'acy of thtT 
method of dsuaUzation. 

A.nother view of changes produced by the splitter plate is found in Figure 
7-l(m) iov'king down on the ground plane, A conventional ground flow/sta^Tiation line 
is formed, 'lut the center of the jet stagnation region Is displaced outward from the 
splitter plate. Again, this is believed to l>e due to the boundary layer formed betw’ccn 
the jet plume and the plate. 

Proceeding then to the closer spacing (S/D = 2) a similar senes of heights 
was examined. At II/D -- 2, a fairly clean jet plume and upwash picture was formed 
(Figure 7-l(n)). At an tl/D - •}, the upwash curved into the jet on the right (Figure 
7-1(0)). A slight mo\cinent of the splitter plaie eaused this to be revei'sed with the 
upwash \ coring to tlie left (Figure 7-l(p)). By careful mo\ement of the splitter plate 
while obsening the oil, it was possible to obtain an upwash flow up the center (Figure, 
7-l(ci)). At laiger Il/D, Figures 7-l(i) and 7-l(s), this sensitivity did not appear to 
exist, liut a stagnation region nutikod the end of the upwash. 

Summnvv of L'lJwasli Pressures 

Flow propeities in the ujiwasli were found by lra\ersing the Kiel and static 
jiiolx'S in the >-divcetion at vaidous heights ( 7 ) above ground. Figvirc 7-2 sliows data 
obtained from a jirohc traverse aeioss the upwash centerline. Note that the static 
piessure was lx;low ambient througiiout the upwash and In-j’ond iis edges. At the edges 
of tlie upwash (smallest value y for v\hleh dynamic pressure is r.ero) and outside Tlie 
edges, tno Kiel pressuie and static pi'cssiue readings vvcic ecpial, in spite of the dif- 
fi'ient shapes and oiientations of the luessurc holes in tlie prolxns. Within the upwash 
tl'e sialic piessuio profile was onlv shghllv deptndent on tlie sire of the static pres- 
s'ute prolx'., The .static pressuix; was nclow amlnent throughout the ..egion iX>twec'n the 
no/.lc c.xii plane and tlm ground i m cpt for a small zone less than 27 mm alxDve ihc 
giouiid lust above the stagnation line. Tins large region of low static pics.sure 


appeal's to be a characteristic of multiple jet ground Impiiigcmcnt flows, and was 
found at all values of II and S tliat >vcre run. 

T);e dynamic pressure profile siiowcd a Gaussian shape except at jet spacings 
below SA3 = 3 where it was distorted by a lack of symmetry. Figure 7-3 show’s data 
oltalned from probe ti-averses across the center of the up\,ash at various heights 
above ground for ll/D = 2 and S/D - 4, For each helglit above ground (z) the profile 
was normalized by the maximum dynamic pressure on centerline (q^), 'Uid the probe 
distance (y) relative to the centerline was normalized by y TMs plot shows that 
data taken at diffei'cnt heights above gi'ound has the same profile shape, which closely 
matches the Gaussian p ofile 

q/q^ = (4) ' 

Profiles of tWs sliapc were found for nozz.e separation distances down to S/D = 3 
providing the nozzle heiglit nbo\c ground was not large enough to provide a significant • 
interaction between the incident jets and the upw-'sh, Ulien the djiiamic pressure pro-- 
file shaj>c Is imlependent of lioight alx>\e gi'ound, at any value of z th.c upwash proper- 
ties on the centerline can bo specilied l»y the values of q J^. and y 

c j • -o 

1 iguro 7-4 shows tlie djiinmic pressure decay along the upwash centerline for 
ll/D - 1 at ciiffeient nozzle spacings. For each spacin’’’ the values of q^ were norn’al- 
ized by tlie coriesiionding maximum uround pressure at the center of the stagnation 
line. Tlie data aic plotted versus ' ‘' ^/ ' > ^ lo account for changi s In the radial flow 
pattern In tlu- ujiwash at different v.dues of S. Tlie data for f>/D ~ and G fall along a 
straiulit line, indicating t'lat Die dvnaniie pressure lanation along the upwash center- 
line can be ie]nesented In a power law decay. As tlie spacing is derreased the data 
de Mates f’com a straight line Ix’caicse of Interfer-'nce bc-twecii the upwash and tlie inci- 
ilent jet How at this nozzle height above grour.d, 

Oui modeling of the fiow dlreeticu in tlie upwash assuivics that the radial 
jiaUeni oi the wall ict flow oii the giouiui continues into the upwasli after the collision 
of w.dl lets ,'t the stagiiation line. Tin; \aiiation of d\nanuc pvcssuio along tin. up- 
wash at constant heigiit above gi ound c in Ik- found from: 
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(5) 



If thei-e Is no intciMction between npwash and inipinf;cment r.onc. Equation (5) ib 
plotted in Figure 7-5 for S/D - 4 and Z/D = 2 for coniparLscn witli data tal;cn at several 
different values of nozzle height above ground. 

The variation of upwash proi>eillcs on the centerline at the nozzle exit plane 
with jet spacing and jet height above ground Is illustrated in Figure 7-G. Decreasing 
nozzle spacing increases the exit plane maximum dynamic pressure, wltli not much 
change in upwash width, until a nuLximuni \alue is leach Ix'twcen 5 and I diameters. 
Fui-iher decreases in nozzle spacing results in a decrease in exit plane dynamic 
pressure. 

At low gi-oiuid heiglits we (.■neountered instabilities In the upwash at a nozzle 
siueiug of i! duunotc'is. v^uch instabilities apjx'ared as abnormally lai-ge nuctuations 
In I'lolse loadings, liul piobc fluctuations In the upwash were generally found to be In 
the lanea' of 10 to 20 p(’reent. I'nstabic upwash conditions Increased the magnitude to 
over aO pel cent. Such conditions were eneountei'cd foi S/U - 3 wlien ll/D was de- 
creased to b'ss than 2-1/2, Corxa'sjionding fluctuations in ground plane pressure \\xre 
noted undei the sanu iinpinuoment conditions, Indliatmg that the entlic upwash flow 
was aftecU'd. ruiilu'r deeivasc in jet spacing provided more stable ujiwash condi- 
tions at low ground heights, but the pressiu'e profiles at the exit plane were somewhat 
uns_\ mmetrie. 

i'he flow pio]-vo’1ies on the upw.’.sh centerline at the noz.'lc exit plane are il- 
lustiati'd in I'iguie 7-7. whicli sJiows ciuwcs of con.‘=tnnt dni.'umc pressure on a plot of 
ll 'n \vi''ii'' .s/n. 'Ihe boiuKhiVR"’ on tne left side of tius cuttc, l.ib-'led 7° and 10° 
lenresent Ihe conditions wtieix' the incident jets o\erlap shear lavors oefore impinging 
on the giound. The lUfloienei' in these two bouiul.mes is the spreading angle assumed 
for tile free let pl'imes, 'Phi^ lower i>oun<lai\ correspond,-^ to a s(.il>iht\ limit iliustia- 
ted on .1 similai plot j'l eseiucd tw ll.ill and Rogers (Refeicnce 8). He tom d tliat stalile 
iiapiugtnu.nl tlov,-' existed it lower values of 11 ^1) th.nu t'u’v indieated, liowevoi. an 
unst.d'lf legion dots .ippear to exist aioimd S/I) - .! as illiist rated In tiu' daslico iigion 
"-kett lied in Figiiic '’-7 
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At the closest nozzle spacinR, S/D = 2, \vc foiuul that the iipsvasli was nilldli’ 
unstable at H/13 = 2-1/2 and 3, but not for lower oi- hl};her values of H/D, Such lnst»>- 
bllitles would disappear if the symnieti'j' of the Iniplngeinrnt flow was disturbed. Mis- 
alignment of tlic ground plane b> a few degrees would smooth out the instability 
descidl^ed above. In addition, excessive blockage by a probe suppotl could stabilize 
this type of fluctuation. 

Effect of Bodies on Upivash 

Since the basic analytic proceiiiues assume isolated up^vash behavior, using 
the tu’o-jets as a building block, a brief experimental study was lUhlei-takcn to determine 
the cftects of adding various simfaces. Figure 7-8 illustrates the effects of the pres- 
ence of our all craft model on flow propoitics In tlie upwash. The total and str.tic 
pressure profiles along the uowash centerline ai-e shown for Il/D = G and S/D - 5 with 
and wltl'.out the model in plact. At tills helgiit aboie ground the change In centerline 
upwash profiles caused b> the model Is significant onlv at points closer than one nozzle 
diameter from the fuselage luidcrsldc. The static pressures increased 
as tl’C flow stagnated on tiie bottom of the model. Figure 7-9 shows tlie effect on the 
upwash of the nuxlols used in the force experiments at tlie exit plane for a lower height 
nboM giound (ll/'D - 3) and a closer nozzle spacing (S/D •(). In t'lis ease the rnoileLs 
altoiod the ilow propeilies of the h\ producing lowered static and total pres- 

sures. The greater the Uxl\ aiea, the givater the etfoet, altlioui-Vti it was not large 
foi an\ ea^es. Tliis effect does nor appear to altei l)n‘ <;eneral Ivliavior of the flow 
but a eorioetlon lor its efteet"- on tlie geneial pressure Kiel mav U' needed. 

r f I e et ol Atlding a Tl.ini Jet 


To determiiie the iinduenee of a third jet on the pixoperties of the upwash 
formed Ivlxu-en two icts, prol>e sune\-- wiu-e tak<-n across one of the upwasii flows 
with and without the thud let in opcratuni. The primarx ut pair was spaced at S/D = 


1. ’J'lu' third ii't was located at eithoi S'll - 3. SG and S/D = a. (iy with lespect to the 


other t"w o lets (3.3 and o. 3) with lesjvet to the line I’cnneiting the first two let centci- 


lliu's (see <5ki'teh on Fignie 7-10 and 7-1 li. 


The impingement geoimdix i.- ilUidnitoil in rigxiie 7- 10, wh.ieh shows mc.isuu'- 
ments for the elostst sj'aeing of tl.i tliini let. I’iie nuVvimum duiamic pres^ule in thi' 
lounfain flow is quite elose to ihe (nedieteil fountain ku-ation. Note that the maximum 



! 


i 


d>^lanIic pressure in tlie upwash formed behveen the prlmaiy Jets is lowered l(f.c l>y 
tlie presence of the tliird jet, even at points that arc far from the fountain region. Tlie 
direct influence of the fountain pi'csencc appears to extend into the upwash about .7oD 
from the predicted location. 

Fij^iu'c 7-11 shows a similar plot of maximum djmamic pressure witli the tliirci 
jet located fuilher from tac primai'y pair. . si.owing that the influence of the third jct 
on the primary upwash almost negligible. The direct fountabi presence again appears 
to extend 1.75D. For this case, probe suiwcys wci'c taken at several heights above 
ground to mvestigatc the fountain flow. Figui’c 7-12 shows a comparison of the 
properties in the fountabi with the properties on tlic centerline of a t\so jet upwash. 

The bicluiation of the locution of maximum fountain dynamic pressure is illustrated b\ 
Figiu’c 7-13. 
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8. JET INDUCED EORCES 


Jet induced forces were oljtained for two-jet configurations by usir^ a strain 
gauge balance for the simulated aircraft surfaces (full aircraft model of Figure 4-1 
and fuselages of Figure 4-3). For the flat plate (Figure 4-10) the force was measurt'd 
withi three separate strain gauge beams. In all cases, only the let-induced forces were 
measured, and the jet thrust was not metric. Jet thrust was calculated from integrated 
total pressi re profiles at the nozzle e.xit plane, including measurements through the 
boundarj' layer. 

Forces on Ali'cratt Mode l 

The aircraft dcsipi used in thi.s study is shown in Figure 4-7, A scries of 
c.spcrimcnts with the basic au plane (Fijturo 8-1) showctl interference force levels 
bocommg moiv negative as the icts were brought closer together, down to a nozzle 
spacing to diameter (S/D) ratio of 2.6, Below 2.8, down to 2.4 (the limit of our 
aiiparaUis) the force ciuwcs stayed essentially constant. The onlv hud of any unusual 
bchai ior was the sUglitly more jxisitivc than e.xpected force levels ami the local posi- 
tu'c tloviation of the force curve for an S/D of 3.2 around a lioight of H/D = 1.0. 
i (Tliis euia e is sliowai a.s a dasiicd Ime in Figure 3- 1 for clarity.) 

St rakes were then attachetl to the lower surface of il.e fuselage These strake 
.surfaces are shown in Figure d-2. >Vhen strakes were added, a much more complicated 
situation de>elopcd. For aid m observing ami understanding the behavior, these data 
arc presented in a senes of three graphs in Figure 8-3. Figure 8-3(al snows that be- 
ginni.Tg at nioderatcU’ wide spaciiigs (S/D =5 7) tlie interference forces ( A I'/^J are 
po.sitne over the range of 11/0 iiucstigatcd. As the jet spacing is decreased, the forces 
become l.irger up to .i spacing of S/D = 1, and then slay apprcvimaiel> coiibiant to S/D 
= 3. 2. I'jguie ti-3(b) (carrying over the curv'^ for S/D = 3. 2) ludic.'ites Unit at a spacing 
of b/i) - 3. 1 the llow bccauu' veiv ur.steadj loi heiglits belwocn apjiroMmately H/D =1.3 
and 2. 1 . 1 lie loiee k\els i.iried uvei Ihe range :in.iicated in the shaded arc.i of die fig- 

iiie. liie bpaeing. deeieab(.'d .sLill lutiher, this unsteady bolia.iior cusappeaied, but a 
it..,inn lit ,ii iiit.^-e '...I re li/F's oriu. leii where t'lc foi ce fell tiom the S/l) -- .i,2 value to 
a lUiich luwiM \ due, i lie width (m !!/l) v:iiii<--') of ihi.s region ol lo'.'-t-.-ed torces grew 
wiui. I a-i .'•/D wab tkeieasi.-d until ior Sjiaeincb less iha.i .^/[' = 2.3 i !'0 ii'-ghei toice 


rcj;ion did not ocx'tir at all. Aftur thib point the i;cneral force level ap^in showed a 
shitht increase \silh closer spaeint;s down to S/l) = '2. 4. Tor these close spacinjis 
another much smaller unsteatly transition area occurred near a height of ll/D = 5. 

This was a transition between the forces with strakes and the lower force levels 
(negative) without strakes. 

The changes in general level and shape of these interference forces with 
strakes were found to lx- a propexiy of the upwash flow Itself, rather tlian its interac- 
tion with the aircraft, as noted it\ Section 7. 

The alic\e c.'qienmcnts were conducted with one specific aircraft di sign. 

L'arlier work iiiiheated that details of the liody geometry could strongly aQeet tlic in- 
terference torces also (Hefercnee (>). An example of this is the effect 1 1 the fuselage 
corner radius. Making the fuselage corner radius sharper produces positive inter- 
ference forces. 

To eonduct a systematic investigation of the effects of body geometry wc used 
a senes of simplified bodies, as well as extens’ve measurements and nioileling of the 
basic .ict,^ipwiish Hows. 

I'w e-Dimcn'-ional lioch rorees 

The lii '-l senes of exinTiments mvolved c\lindrical bcxUc.‘ (.shown in Figure 
I-;'!, siiiular lotlie aircralt fusi. lage in flu' previous scctiui, with var\lng lower sur- 
f.iee eontour. In .iddition to eonliimiag lh.it a fuselage-lvpc surface can experience large 
positi\e toices in-gi ouud-effeet it w .is found that there was a eli.mge between .ittachcd 
How .iiouiul the Ixxh and so]\ii.ition .it the lower comer as the comer radius was 
ch.mged. Flow \isii.ih/.ition pliotos us uvg tufts show this Ixlun iov cle.nl\ (Figure 8-1). 

'1 i'c .s\wti. Ii beti' een att.iehed .md st>p,n .ited How w.is four.il at a corner radiiLS to lody 
wuith [\/\\) of tk ISS. -\t this i.idius both flows could c.xist, .md could i>c switched by 
distin I'.nu’f's of the flow. 

I lie eltect of comer i ulius on llie intoiferenco loiccs is shown in Figiire 8 - 3 . 
I'HUirrim lioin a sii.ij.]! tointi winch icsult.s in .i l.n-gc pos;ti\e intcMci cnci’ force, 

Uie foice dec re.isfs .is the eoiiici i.uhus is incrca.scd. When the flow switches from 
SCO 11 itioii It tl'.e comei.s to .Ut.ii lied f1c>u ■' l.sSl, tlie loiee lumps iipwaid *o .i 

huger \alue. 1 uiliu t incrc.i'-c in the radius losults in a ictuin to tlie tiend of dt— 
cie.isipu lone with iiii le.ismg i.idius. 


The occurrence of a larger upwa’xl force (liody "drag" in the upwash flow) for 
attached flou than for separated flow is understandable with reference to the pressure 
distributions of Figure S-G. The pressures on the lower surface anil the upper surface 
aie affected by the change from sepai-ated to attached flow. The pressure on the upper 
surface is very' near ambient when tlie flow is sepai-atcd. When the flow is attached 
the ujw’ash above th« body acts as an ejector and produces a larger negative pressure 
which, in this particular case, is a positive interference, or lift. The prcssui’c 
distributions of Figure 8-6 also show tlic reason for the sensitivity of the forces to the 
body geometry. Tlie integrated pressures (foi'ces) involve two areas on the bottom 
surface, one positive, one negatne. The resulting force is thcrxjfore the small differ- 
ence Irotwcen two large numbers, and very sensitive to shifts in the curve. 

Another study area was the effect of the strake depth on interference forces. 
Figure 8-7 shows a continui'ig increase in upwards forces with increased dcpUi. The 
curve for r - 0, no stinkes is a rcfei'once line. This was the curve of greatest inter- 
ference force in Figiu'C 8-3. Nolo Unt also a boi’y with acorier rai'.ius and fittec) with 
strakes wluch arc level vvltii the fuselage l»ttom section produces a higher interfer- 
ence force tlian the flat bottomed bodv. 

With a stralie dopth. of = 0.125 a series of exiieriments was conducted 
with throe difieicnt Ixxiv widths, U/IJ = 0.5, 1. 0, and 1.5. Results of this senes are 
shown in Figure S-S. A piogi'cssion of higher loi-ces with greater body widUi b seen, 
bia the thinnest IkkIv does not have as fast a drop off of force witli height alxjve ground 
as found for the other two. All thiee of these txMlies are enveloped by tiie upwash for 
values ot ll/D down to 1. 

Rxtendmu this investigation to other classes of Ixxiics; a lower fineness 
ratio (L/D - 2, 8 vice IG for tlie fuselage type Ixxiy) rectangular Ixxly was tested with 
and without straKcs (Figm-c i-8). Nozzle spacings uetwocn 2.4 arid 1.0 result in nozzla 
contained within tlie planfomi. 

Intel fereiKc toices on tlicse Ixidies vvitii strakes (Figure 8-9) have much 
l.Ugei negative values tivan tliose tortile in cceding cases. Tins lesults Ijceause of tiie 
laigc 1 ;.iea lor the lower si.u f;u-e pressures to affett. Wo found no signifu .ml ctfei-t 
of t'lc eonier t aduis or the upper surface corner geometry toi tiiis Ixidj. 'i he addition 
of stiahcs agaia pixiduccd a largo jiositivc increment in foiei'S, and a region of 




unsteady flow de\ aloped for S/D =3.0 and height above gi-oaiid between S/D of 1. 5 and 
2.5 (Figure 8-9), 

A survey of tlie effects of strake heigirt j nd s+rake si)acing were conducted at 
a nozzle spacing (S/D) of 4.0 (Figui'cs 8-10, 8-11). Variation of the strake spacing pro- 
duced some changes In the force curves, with an optimum at some intermediate 
spacing. In all ca/ as significant jxisitlve forces resulted. The positive force incre- 
ment due to the slrakes (Figure 8-10) Increaseti with increasing strake height up to 
alxmt ^/D = . 25. Doubling of that \:ilue produced vcr>' little change. 

Also examined briefly was the flow with a ver>’ bugc plate at the nozzle exit 
plane, using a square plate 12D x 12D, with a nozzle spacing of S/D ~ 4. 0, Lai'ge neg- 
ati\e forces were foiuid without sti'akes, but for this 'h'chlcle” the large positi\e force 
incremerts with strakes did not occur (Fijrure 8-12). In fact, for some values of liAl 
the strakes resulted in more iiegatue forces. 



Sn m*)ol 



1 I fi I ( III u't ''HafiiTj, (VI .uu'i'ifi force' 


O + 













IU.-C '-I. Hoa, \i9iu\li/uUon of up\\.i<?h on fuscl.i'^c. 


N'0-!/l.CS ^ 






*T53S»>Wl 


Pv 




B-mW JdiHItf' L/0 • ftO 
J«l S'O • 4 0 

Strike if/f> • 0 !?5 


Boify vwdih^'Svf^'tH'Vt 












9. THE DE\'T: IDl’MENT OF A TIIEOHETICA 1, MOEE L 


The problem of single vertical axially symmetric jet lmnln>;ement has been 
studied Ixjth cx]3erl mentally and tlieoi’etlcallv by many ln\ cstlfi.itors. A survey of 
some of the pertinent literature can be found in Ref. 9 witli a sketch of the basic Row 
shouai !n Fjr. 9-1. Do\^•nstl’cnm of the nozzle there Is a velocity decay region due to 
turbulent mixing followed by a stapination region on the ground plajie and deneetJon 
parallel to the ground plane until the pressure returns to an^bient. At tliat point the 
velocity is a maximum. Turbulent mixing will then goveni the bcliavlor of the raillal 
wall jet. 

The impingement region has lx;en treated theoretically for an Invlseld irro- 
tational or unifonn mcident profile (e. g. , Hef. 10, 11, 12). An amal>lic inviscid solu- 
tion (i.c. , Euler Eqs.) was gcnoratcil for tlie impbigement region due to a rotational or 
full\ de\elopi'd incident piolllo in Ref. 19. More recently, numerical solutions liavc 
been generated, for tiie lotational uuisciu impingement problem using a frozen vortlc- 
it\ concept (Hef. IR. These approaches, although theoreticalh sound, have dcficicn- 
nes wlien comtuircd to exqwnmental d.da c\en when coupled with a boimdary layer 
.malvsib, I'oi’ e\.iirj)le, .'n\ imiseid solution of tlie imjnngemont region '\ill yield ; 
w.iil jel giowil) tb.il is inaersel> jiropoilional to the i.adial dist.mee fi-om the stagnation 
fiouit .uid .t niaMinuin \eloeiU along the wall equal to the imixinuim wloelty of the lu- 
i-ideut profile. This U;h.i\ior is uo\ei quite aclue\ed in the physical llow. No ade • 
qu.ite theoietie.il .m.iUsis is as.illable that couples the jet deea\ , impingement region, 
am! result:'.. it turbulent w.ill jet. 

In the present stud\ semi-empliie;ii ino<lels .ire fonr.ubiti'd for e:K-h iv>rlon 
:u',d eoigileil togetliei to snnulatc the plusu-.d tlcw. Momentum Is conseixed .iiid the 
empirieal input to each legion Is well defined .iiul eapabK- of cxqwnmunt:!! \ enrieation. 

'1 he est'ililishmeiit of .i \ ertie.il let Impingement model is the basis toi Ihe *k'- 
luiMoi ot tMO lets imp.nguu’' on .t gronmi plane. A sketch ot the Ooa pioblem is shown 
in 1 igiue ‘1-2, riuiue u-2{.0 •minors the \eiUeal o! me .■.ip.iaiemg the jet siagn.aUon points 
.':ui the givun.l pi me. )i l lmpi'ig<m o.i the giounl pi. me and detlfct'- to form u.dl 

jets. 1 '>e w.ill jets llun iinei;iei .ind form ,m upw.mh i!elli.‘etion /one where ‘lie w.ilI 
jets Lollid.e ■m.ii lie tui’ii.(l I'pw.mi le.uiii,; the gioumJ jii.me. Ihe st.ipi.ilion line lie.s 
on the ground in the M-ilieal pi. me of simmetiv •let'ween the two niijuiigiiig jets. The 

>.-l 




& 

f 

niiL'anium mwabh stupiation pressure occurs at the midpoint of t!ie line connecting: • ^ 

the tM'o jtt stagnation points on tlie i;round. The ttrounci pressuie tlicn di-ops off with 
latci-al distance aloiiK tlie stagnation line. Figure 9-20>) shows tlic radial stre.imllne 
pattern that has been observed (Kef. 3) both in the ground plane and in tlie verti- 
cal plane of symmetry between the two jets. This flow situation only exists when 
the jots are spaced far enough apart so that the jet impingement zone does not 
interact or has a negligible effect on the upwash deflection zone. 

Figure 9-3 slious a sketch of the flow pattern when the jets )u.\c a spacing 
siioli tliat the jot Inipuigoineiit pressure rccovcr> is not ai'hietcd Ia'idi-c the ujws.isli 
tli'lloition occurs. In tins situation, some of the wall jet streamlines interact \\ith the 
upw. isli liefleetion /one so that ambient pressure is not aelileteii along the entiu- iiji- 
wash deneition hue, lleiKO, a iKiundan of minimum pressure occurs bobsetn the 
jet Iminngement and ujiuash deflection zones (Figure 9-3j. Tlas boundary is a resell 
of the matching ol the prt'ssures in these (i\o legions. 

In general, t!ie flow in the upwash sheet l>e>u\os as a free shear flow in ihat 
It cvliibits t'-i'ieal tuiiiuleiit miMiig jeluuior tor tlie decay of total pressure and \cloc- 
ih , I’lgiue 9-1 shows a skctili based uiwn emjiidcal obsen .itinns ui what ma\ h.igj'on 
wl'i'U Uk' two jets ai-e ><■ r i elo.ic oi uisl jmor to eoalcsi'inicc into a memed itt witii ro 
iipwa.-li thm , Lpw ish How dots not exist in the \eriical plane containing Ihi i-'t stag- 
nation i'omls I'ut doe., oeeu! l.'teralh, 'J'.ne sticnmlines are (iiieiterl alxjut :i rone of 
no upwash llow near the eenlei ol tiie iij'wasii sheet. 

It is the intent of tliis .stiuK to mo<lel the ix‘lia\ior of two closel\ siiaeed jets 
into I tine t>.' lo; in .ui ni>wasti sheet. Few imi'stig.-.toi's lia\e attemjiteci to stud\ this 
pivL'leni ill ain det.ul (Hef, 3, 1, J(> and 17) and little was known ai/out the Ixahavlor 

ol i'l'"~ei\ f-paeei! jets jn ior fo tins studv. 

Aliliouga ^oine lesults for f'-.e tw-o jet mnnnueincnt proldem have luwn nb- 
lainoil l>\ auiiK’; k’al iJileualion of tlie Naviei -Sli7ses eriisitioi'.'' . (Itef. 171, tlii.- !\po i 
ipriuaeli I eitauii\ i,c: .t pir.ctieil .(pjnoaeh at tlie turieiH la^'ie. Turbulci'ee incv'el- 
lUu, ^ I'.an in leiu ih .'e.'lt " . am! bound n » soii'luioe-., u-^uii fi ruu t omjH.a.iti' la.i! time, 

,'...se ' li >; iMiPiliie t .',r III ol it "u 1 11 i.g vallil loll', o: ibu ‘'..'iVit!' •> • ii,.i 

iioii'. \ ' .■ae'i s '.luiii -einpii’ical is>ai,'>‘ u.*, iit imuii'liiii a,'.'iionn t iKe'< 


Jet Liipingcmcnt Models 


The j* t Impingement problem is subdt\-lded into three major reglop.s: 

® Free jet decay 
e Impingement or deflection 
• Wall jet 


Each of those finite regions aif farther sulxHvided into minor regions where 
the flow transitions from in\isoid to turbulent behavior. For simplicity, a power 
function of the form, 

fO;! - (I - 


is used to simulate lx)tli pressuie and \clocity beha\ior. This functional form was 
fli'bt used by Abramsvieli (Hof, 18) lor shear la.NX'rs. 

h. 1 Free .let Doeav Model 

Figuio 0-5 snows a sketch of the Ivisie jet decay mcHicl, suUil\ided into three 

rcjiions: 

• potential core 
o tiaiiNition 
o fidh tiu'l'ulcnt 


The half \eioeih> boundaiy giowtli In each region ean Ix' determined lai-geh 
fiom exiierimental dot... T l.e nuxiel assumes a simple linear I'oiindary !^l■owlh in each 
1 cgion. 
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The constants i^ovemlns tlio half velocity boundai-j’ growth In the transition 
and fully developed I'cgions have Licen determined cminrically. In addition, ihc 
potential core length and fully devcloiMid similarity exponent of the velocity profile 
must 1x2 specified. 


The prowth constant of the half-velocity boundary in the potential core 
rcfplonis determined by matciiinc: the boundaries at the end of the potential c*orc. 
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In addition, a linear decay In potential core radius, r^ , is assumed, or 
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ruilher, the paniircter Is determined by nn telling tl\e Ixiunfiaries 

tM2tween the tiansition and fulU de\ eloped regions. 


'I'l'.o static pvcssuic in the jet ai tiic no//le exit i:- as'-umed to Ix' amid''nt. 
TIk inoinentum equation ran t. hen l>e n utten as, ' ' 




-n I I 
■ u ' 


J diT 


' rdi dv' -• t'0,\b.'l AM 


( 9 - 1 ) 


let \ (i ■ /» 

\ r/) 

III 




'1 . 




! ')-! then reduce^ to: 


V, he 1C 


/ ^ \ 

/ "i \ 

-- 1 

- 

\ \ y 

\ • / 



1 

1 

> 

1 

\ ) 

1 - .1 

Ti.) 1 


S>\ 


1 


( 9 - 2 ) 

(0-d'i 


n-i 


Cj.., Is nn Intcfji'al of the profile function. In the fully developed rct^ion 
J 

this Integral is assumed to Ije a constant. 

The general form of the velocity profile is chosen as, (after Ref. 18), 

Of, 


-L r , ( . / ”■< ' "jr. \ ■’ I 


(9-4) 


where 


'JC b. 


Substituting into eq. 9-3 and integrating, yields 


>) 


'bv/’Mc' 


f 1 ? 


",r- 


0,1-1 3o , 1 2 lo . ) i; ^ 


1 -i— i 




_j _ i> i ^ 

) ■ 1' , ♦ 1 ' 2a *- 1 " 30 ' 1 -lo f 1 i 

* J li J ' 


(9-5) 


whu’h IS a general OKiii'ession for 0^,^ y integral as a function of potential core 
ladiiK-' and \’clocit 3 c'-’jonent o j. 

The ratio of iialf \eiocity widtli f jet width can be determined from (9-1) as, 
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9. I. 1 I’l-ileiilKil Cmc and 'J iMiisUmn Region - Tlie momentum equation 9-2, given that 
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Thus, given v and the half velocity lx)undai-y bjjj, the solution of equation 
(7) yields the exponcntoj go\eming the \elocity profile in the potential core region. 

The exponent Ojp^ is known at Z = Z^, Tlic fully developed exponent 
Ojpp is also ppven empirically for Z > Zpj^. 

Thus, the velocitj exponent mustvai-y accordingly In this region. By 
prescribing the condition that the jet half I'adius (b.) must vary continuously from 
the end of the potential core x'egion to its known \alue in the fully developed region, 
the ti'ansltion exponent n-j-p can lie computed as 
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Thus, the momention equation (9-7), given n jq ~ 0, can lx: used to determine 
the maxi mum velociU : 
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9.1, 2 I'lilh Dexeloped Hegion - In this icgion similaiity of the )( t profile Is 
assumed ( jj.-p ~ 1. 1>) and, thus, the ratio of half \clocity width to jet radi us 
and the siimlaiiU integral ate constants. The \eloeity, as determined from 
equation (9-9), decays imci sely propoitional to the half width or 

I'lCLire 9-6 shows a t>pical I'esult for nreximum velocity decay and 
l)Oim(lar\ giowth in tlie fiee jei compared with some test data from Hcf. 19. 

9.2 Jt'i Dc'lettion Uenion 

Tiguit 9-7 show"- a sketch of the seating jiarameters foi' the jet imixiniiement 
problem. At ■<ome height, Aj, the incident let "feels" the eifect of tne ground jilane 
ami the einteiline \eloeil\ decieases and linall, slagimtes at t''.e go omul. The jet 
mixing IS inmliited sonn w i’.'t in the dei elei .<t:on lepi'in, althougli t' e stegi'ation 
piessure is olisen. eel te) le !ov e: then th.e total jnessurt of the jet tt the ground 
elUat in'iglit A due lo tue eificis ol iii!\ing in tile lieceierat'or le'gmn. '1 Oe 
impnigenumt iialt xcIocUn i adum ami halt v. idth w ill I;, usee! to luil’ncr s< ale the 
wall let legions. 


The groiLid effect height has been determined experimentally and is gi\en 
from reference 11 as, 
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The stagnation press lu-e is also given in reference 11 as, 
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Tlio id rnpingemciU raciu arc determined fi*om the jet decay model. The 
iniinngemenl or deflection radius on tlic ground, r^, is defined to bo the location 
whore the jiressure in tlic ground jil.ine reaches ambient conditions. 


9.2.1 riroiind Pressure Distrilmtion. - The empirical input data rcciuirements are. 

' Ilir 'o’ ^ ground pressure recovery function is assumed to hoof the form. 
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A momentum balance is performed with the integral of the ground press ui'e 
equal to the momentum or thrust of the incident jet. 
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This equation \ iclds tlic exponent <tt . of the pressure i'ccovcr\ fimction that 
satisfies tlie momentum condition. 


The deflection r.idias is estimated empifically and in nonciimcnsional 
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Figure 0-5 shows the kind of con.i'^ion achieved for the pres'^uro 
distiibiiLion on tl'.c ground from this sinijile proc> ' ’re. I iifoiliinatch , the degiec 
01 ‘'Ucetss larc'-h ileponds upon the cliosen form oi the pressuix; uiiution, 
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9 . 2.2 Wall Jet Transition Model. - The \sall jet transition ino<JLl coreJsts of thi'oe 
subregions: 

i) Inviscld deflection region where the effects of viscosity are assumed to 
be negligible except in a region close to the vail and near the edge of the 
deflected flow. In this region, the inaier boundary layer is established as 
goven’icd by rwially symmetric stagnation flow. 

li) Transition region where the effects of \iscosity are Ijcginmng to 
dominate and the inner boundary layer and outer shear flow transitions to 
the fully de’"eloped turbulent vail jet. 

ill) Fully dewloped flov' where the effects of turbulent viscosity dominate 
and the static piessure tlirough the wall la^er is eonsidez'cd arnlnent, or 
fully rcco\ercd. In this region, the nearly similar wall jet dc\elops. 

The inviscid deflection region is assumed to occur u der the half veloci.y 
vidth (I'yjjj) of the incident jet proOle. The preoS'irc distubutlor on the g'round was 
determined in a prcMOus section as 

o 

AP -- ^ t - (7,^^) -J ‘ 

The maximum \clocit\ in this region can ize predicted adtquatt. 1> from 
IJcrnoulh's equation as, 

' M , '' j'|. 

— v.here 0 ^ ^ ~ (9-12) 
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Tlie initial Ix;undai\ la\ cr thickness i ^ at the stagnation point is assumed 
to IX- ix)th a fimrtion of nozzle height alx3\e giound and of Re>nolcls numlx-r. (Ref. 20). 
6 ^ is .issumed to Oe ilctorpiipeHl CMi|)incally . the detcinin.ition of Uie "laxiiiiuni 
cil'xitv vanalujii i'- the oiiiv requ 1 reini lU of this region since the .■na\iinuih total 
pi i' ‘'Til i c 1 '- assiimeii coiist;inl. 

\ nrxiel must 1*. L^tabhsnt-ii tliat rn.itclu's the liu i--.ci'l m.i .iinuin -•olocity it 
toe end ot the must id legion and 1 1 .iiisu iops to .1 lulli, unhulent .di |et. I igurc U-9 
vp.)-'.^ ,1 sketeh ol tlic wall jet tiansiUon Izoin'daiv gro".!!' n''xi< 1. 'i I’e niodol cou-'ists 
of tvo l.ieeis. 'iicc inpti ixmniiaic lavtrllov. and outer slu ar. i'e fiil,’. di- eloped 
!>. ha\ ie>r of tne IxniPd u v la\e r and i'al f \ eloc in w idth !>■■ .e Mined to • f , 
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For simplicity, Pg and 1« 

The form of the boundary layer prowtli arises from ex]x,rini. ntal observation 
(Hcf. 120). The tluckness of the boundary layer appears to relatively constant in 
the vicinity of the stap^ation point or in the inviscid deflection repcion. 'I’he ixnmdai-y 
layer prows relatliely slowly until pressure recotcry takes place and then the 
iioiuuiaiy layer pi’ows with radial distance from the stagnation point. This delay in 
boundary later grwvfh is tlie reason for the choice of the aixite function which 
precludes the establishment ot complete sinnianty in the wall jet profile. 

The empirical co’tstants for the fully dcteloped wa.'l layer gcravth arc 
approvini.itcU , (lief. 20) 
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Th'is, fii ' limit o'’ tl'.c r.atio of tl.e lx»undar> layei 1‘Oiglu to the haif velocity 
!>oint in tl;c snoar l.nor :s, 
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Tin's '..liuc IS in agreement with other investigators (Kef, 21). 

I'bo inofe .im’.mroons value of rif, in ti'c limit can i;c computed from. 
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wl»ich is sm;ilU'r tlmn l!io often quoted v;ilue ol i/O find mny be the 

result of the tinito tonii of the ehosen profile. 


9.2.3 Veloeltv :tiui Stiitle Vi e.ssure IHoflles , - I'itii'.re 9-10 shows a sketch of Uie 
two layer motiel chosen to rejneaent the velocity and static inessure profiles In the 
dcveloplni; wail layer. The inner Ixuindarv layer velocity profile is appro-vlnutcd 
by the well known turlmlent power law. 
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The static pressuiv across the houmlarv layer is assumed eenstant, or 
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i'he i.tlui of the halt \eloeitv lua^ht to the wall laver heieht I'.m lie 
delenniiied h\ , ' 
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‘I h. ! U uiial Moiiuntiiin Mu\ l.ijiiition. - The aowiniu;; eouauea. toi tlie flow in the 
all III i-' lormui.iieil s', i i.tdi,,! mnnientum Hiix equ.itiiin 

1 ’.' aa !ii\i''iui How seUi*um o! tlie e\a».t ii-t itiinnd eqv.uu'iis (Itei. 1 1 ) il e m 
U- --he\ui, itnin the t:n UeM .■aial\‘.ie,.>l •.oluliun, ll- ti ihe maiauiuoe ol tiie tiiooientun 
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flux in :iny \v:ill n‘t sector summed o\cr 2 n is exactly ec|ual to the thrust of the inci- 
dent jet. Tills is the t.uslc prlnclnlc used *o c-onstruct the folloixlng equation wlicrc 
a pressure term his Ixien added to account for the static jircssure variation in the 
deflection I'eglon. 
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The :ilio\o cfiuation applies lx>th to the transition nqxion and fully developed 
wall jet. In tlu' traiisillou re>;ion, a static pressure \anation exists in the wall layer 
due to the streamline dcneetion and hence the pressure term on th.e RllS of tills 
equation, llecaiisc of tie stioamlinc deflection, the momentum Hux in the ladial 
diri'i'tion iloe.s not ai'hie\o its coiishint op.tlnnim \alue until pressure recoveiy occurs. 
The function LUS reflects the change of radial momentum in the 

deflection I'cgion. 

Hence, the l-cha\ior of this function must !x? suth tliat, 
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Tiie si'cond luiictum f^ on tlie I.liS of this equation is a term lenccting 

•■'vis ' ' 

the hisses in moiiientmii due to skin iTicUon or shear at ti'c wall. This term Ix'comes 
ii ipoi-t..iit ai lame distances trom Hk- stagr.ition point. Its behavior c:ui Ix' .i)'|iroxi- 
mated in , 
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In p.ii .iiiU't rie tviriii tin.' radial ■noMctiUiiii lla\ eq . (9-17) can he iurtlier 
I \)'i e-'-ed as , 
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The integrals of the velocity squared and pressure profile functions vary 
with radial distance, or 
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Integi'atlon of the profile functions, 9-11 and 9-15 yield, 
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In tl’.eso equations a ami K. are functions of r. 
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The nminentum flux equation can then ix* rewritten as, 
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A'' inentioncfi mvMOusIy, the ladial niop.cntuiri function ) is relate 1 to 

the nies-5iiie' i« i-ove:} m the dene'etion I'egion. The approximation is r.ow intiodueed 

tnat tlu ladi il momeetum llax is lelated to the )>ies-'ure reeo\ei\ tunetion I ' in. 
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The behavior of this function Is sucli that u heti pressure recoveiy occin-s at 
the wall all of the streamlines ir the layer arc assumed parallel to the wall. 

The maximum velocity at the start of the wall jet at r = rjj|. is (j^iveu by the 
Inviscld nernoulli equation, 




(I - f (n.)> wnere 

tr « 




.IHi 




The Initial lioiindaiy layer thickness 6 ^ ia assumed constant throujjl’.out U'c 
linlscid rcj;ion, Tlie initial height of the wall laicr is estimated to be the gi'Oiind 
effect height A j« Thus, has an initial value given l>y. 
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Gi\en tlie inuscul starting m.ocimum \eloeit>, tlie momentum flux equation 
can K> sohed for the initial half clocity width I'le pixifilo as ' 
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FlRure 9-11 shows a compai-ison lietwcen the Initial profile predicted by this 
method and the exact solution for an inviscld ix)tatlcnal fully de\'eloped ImpiniilnR jet 
at SOUK) distance beyond the deflection region as given In Ilcf. 11. Tins simple pro- 
cedure yields an inviscld starting profile not unlike that of the exact Inviscld solution. 

To compute thi.* maximum velocity variation and oouiidary along the wall, a 
fully developed value of the exponent c must tjc prescribed. This value is chosen 
to be = i. 5 after Abramovich (Ref. 18). The turbulent wail exponent N must 
also lie prescrllxid tlu'ougliout the wall later and has ejcpcrlmenUiI values ranging from 

7 < N < 15 


'I'he \alue of 7 is chosen to lx; the lully deteloped \aluc at the end of the Impingement 
legion. Tile tadation of the e.vjionent througtiout tlic deflection region has to lie 
prescnlx;d wiUi the boundary values gixen l>v, 
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Due 111 .1 latk of cnijM ru .il or theoiclical etidcnce in regaixi to the Ix'iiatior of the >uo- 
file in t’lf deflection lei.ion, ;i simple hnea’' \an.mon of the exiionent \\a.<5 as'^umed, 

Tiic li.slf \clocit\ lieight is deteimined for the initi.il profile .md its U*ha\ior 
.it tlie end of the deflect. on region is al.so known. Hence, the half \clocity ix;h.i\ loi 
a.ssumes a simple linear x.iii.ition matelu'd to tlio full, dexcloped height at Uw end of 
tlu‘ (iefiection region. 

'1 lie m.vxinuim \elocit\ in the dcllcction region c.in then lie computed as. 
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In tiie fully developed wall layei*, this cxin’cssion reduees to, 
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Fiy.uic 0-12 slious a typical numerical result for tlie boundary layer thick- 
ness!, as a ratio wltii rcsjieet to Uith tlie wall half-velocity tluckncss and layer lieljj:!;t. 
I'luurc 9-i:i shows a twplcal \ aviation of the wall |et profile as a function of radial dis- 
tance. Tlie \ aii.ition is ie!ati\cl> insifpiifleant and becomes lery neai-ly similar as 
the radial disbinec lucre. ines. I'iKUire 0-1 1 shows a t\^ 1 lcal oierall result from the 
wall jet tr.msition model in terms of la\er thu'knes.scs, ixround prcssiu-e, and nuLvl- 
mum \eloi.it\ .ilou;; the wall, riy.ure 9-15 shows a comparison with the data of ivf, 20 
of tJu' m.LMimun lelociti \arl.ition for different no/,/Ic liciifiits. Uenoralh good agiee- 
iiu'iit i'' .icl’.ciied. In comaanson, the lc\el of m.ixiimim iiiviscid ground velocity 
.it each no/.'lo luighl i.s ,dso m.ui<ed oil. 

‘ 1 . ,i 'l^^ o-.Ie t In tel m tnm Modi 1 

I ignic 9 - id • ho.'.s a sketcii ol tiu ground jilano cixjrdin.ite systems lor tlie 
P.u) H't iinpinj'eineiit inU raction proldem. The lets sf..gnate on the ground, cii-flect in- 
to w.dl jots ,iiu! I'olhdc .dong Iheir jd.iiK' o! s\minetry. The lorm.itum ol mi 

upw.i'.h si.igii.ition .'lino i> .i result ol the upw.ird ileilei’Uoii ol tile w.dl jel ihuv. 
riit' point Ml. on iho line n'lning the ml oi'mei'N. is the onl\ triu' stagnation 
poinl u!'i’io the Ilow vsmn'^ i.) lost and diMloets \ ert le.dly . \t juunts oil oi t!'e 
ooMt'M . iMi' flow iv liolieeos: .It .III .inulo such th.it the pri'ssure is i-epresentativo 
oi the '-l.iui'atuni .M ihe norin.d i-<>'nj>onent ot the iii.iMinuni ill let \’oU>cii\. '1 he 
Ilow isM'tinui's r.idialh alU’r d-hleetum into the upw.ish --hei't is li.is been ob- 
•-o! \ oil o\poi nn.'n'alb 

'b i 1 M iM uvUii i’le-^'Oio in'-', tabu lion .dong 'he I'pw.i'.li ri.u.isilion lino. lo 
o''luiu‘ti ill. upwi-.il I 1 ^, n ip.'i I 1 ' u’ 1 ! Oil'!-! ;ii o-.--ui o d. '.I j i.hiilion 1 ''o m >; la.ii 
\oisu\ at •no 'l"p’U;.o! hPi ainsi l>o ...Jiapuioil 
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The ;irfoSuir; is then, simply. 
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^^^ure 9-17 shows son\c typical icsuits comparccl witli test dntn. for tlio 
maximum .sta)^natlon point pressure on the j^ixiund (as computed from the wall ict 
transition model) as a function of jet spacing, Kairl\ i^ood ap'cemcn^ !s aclueted 
down to a )ct spacing; of two diameters as a result of the transition rnotlcl. For 
comparison, the full. de\ eloped wall jet decay is evtrajiol.ited hack and yields much 
iuglier pressures tlian t!ie test data indicate. 

In ucneial, the wall let transition motiel must he used to eomjurte tlie 
maximum \eloelt^ alonp tlie staunatioa line lor close jet sixuinK. For a jet spatinjt, 
s 'll, lai^e enoi 4 ;li such that tlie wall jet is fully de\elojicd .md the |et deileetion 
regions do not interact with the upwash deflection /one, a simple exjircssion can U- 
demed to )u edict tlic iraxinium jnessures alonp. the staijnaiion line. 

Under tliese conditions, the wall let m Lxinuim \elocitv can Lhi wiitten as, 
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If this expression is fuilhcr nonclimenbionahzcd by the stagnnuon point 
pressure at x = 0, ,^p 
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oi, since S^ =~ S , 
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rmiiie 9-ls shoMs the computed results in comparison witli test data for 
si\cr.il )et spacinps. ’I'lie alien c e. 'ivssions cor-espond to the lari^cr jet siiacir.f^s, 
1 01 the closer ict sp.icinj;s , cciuation (9-27) is inodiried sh<,ditU’ to account for the 
■-lalie piC'-siue laiiation in ihe deOeetion lemon, or. 
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• 01 , ’M\-> 1, ilu alio'i txpits'-ion lediucs to ti',e_ onu,inal cciuation (',)-27'. 
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The correlation Is very good for spacings above 3. 5. For closer jet 
spacings, the model yields slightly higher pressures along the stagnation line. 

0.3.2 Model of U pwash Deflection Zone Without Jet Deflection Zone Interaction. - 
In this section, a model for the ground pressure distriliution of the upwash deflection - 
zone is developed for the problem of lai'ge jet spacings. WOien the jet spacing is 
large enough, such tint the jets deflect into radial wall jets prior to intei-action, an 
upwash deflection zone of high pressures develop. This zone is ricmarcered by n line 
of ambient pressures. E.vternal to tliis line ambient pressures exist in the wall jet 
and Intenial to this line (y<y ) the wall jet bcfpns to deflect upward and, hence, 
de\ clops the upwash deflection zone. 

Ti.i* gi'ound pressures in the upwash deflection zone exert a force on the 
ground plane equal to the normal momentum of the entire upwash sheet. The noi'mal 
moniciitum of the uiiwash sheet, assuming a radial streamline pattern emanating from 
the jet stagnation point, can l>e computed simply as 
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foi cue ([U.idrant oi the imwasli sheet. In this expression, Uic momentum oi the inci- 
dent jet is ilnided e eiilvjKr unit radi.ui into the wall jet and fin.illy into liie upw.ish 
sheet. The normal moinentuia of the entire upwash sheet is then equal to 2M , 
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oi, roiighli , , of the momentum ol one of the jets, 

I qiiatiug the noriv il momentum of the upwash sheet to the upwash zone pressure 
(list 1 il'.rtion, Meldsthe follow mg expiession: 
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Relating thib expression to tlic cylimlncal ground coordinate system of an 
incident jet (see fig, 9-iG), and noting that 

\ = S , tan<;f> 

2 11 

2 

iL\ = S , scc~(^ ^\(^) 
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the above expression incomes, 


or, 
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To de\elop this ex-jiression linlher, a pi-essure function must be assumed 
aeros"' the upwash 7one at any location or angle0u. 
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'I !ie cxines-Moii lui A the stagnation line v.as de\clop<jd prewouslv, 

i|, 9-Js, Milistiluting all this into eq. (9-.}a), and noting that 
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dlbtnbution. 

Now, 
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let C = f(n )dn be the profile integi-al for the gi-ound pressure 
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An hypothesis iibout the bchiivior of the upwjish zone must now be intro 
duoed into oq. (9-38), as, 
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Tins cvtiressioii says that the upw.ish zone becomes iniiiutc in width and 
the radial distance to the stagnation Imc also boeomes infinite as 
Hence, the mass in Uic wall jct also approaches infinity. At the same 
time, tlie .stagnation line pressure approaches zero. Substituting cq, 

(9-39) into tlic integral expression cq. (9-3S) .yields a finite intcgi'al. 
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To deri\e a ^im; 1 le expression toi the upw.isli zone, using eq. (9-10, and 
m.ikmg the fuil'iei .issumjition that there .ire no momentiun losses in the wall jet, or 
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ihe lollowing expr'.s'-ion t m ne iL-^ed to ostiP'-ate the width of the iqw\ush zone. 




.1 

sO, \ \\ J 


(0-11) 


u 1 , , V K . v' 

\\ I 1 1 -'ll.; U >11'. 


Thus, the initial upwash width is proportional to the jet spacing and depends 
on the similarity constants for the velocity squared wall jet profile and upwash gi ound 
pressure profile as well as the slope of the wall jet half velocity growth. If it is 
further assumed that the exponents ^ of the ground pressure profile is 1.5, an 
estimate can be made of the initial width of the upwash zone. Assuming a fully de- 
veloped and nearly similar wail jet flow, or K , - . 0962, = 1,5, and 
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For c\am))le, at a spacing Vd,^, = G, the entire width of the upwash deflec- 
tion zone along the lino connecting the jet stagnation points is appixiximately 2.4 
nozzle diameters, A more realistic or rcpi’csentative boundai-y for the upwash de- 
flection /one might be the half pressure boundarj which can l>e easily computed as, 
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9.3.3 Ipw. sh ?iIoiacntum Models for Close - Jet Interaction. - To fuither estimate 
tiio l>cha\ior of the upwash zone on the ground and ujiwash sheet for closer jet spacings, 
.. mexk 1 foi t'.k. iMzliacior of tl.v momentum in the uiiwash sheet must be inpothcsiz.cd. 

U'.e ii.U" arc >i)<tCod clo'-er. their boundariCi> bognii to merge, and a redislrihii- 
l I'll ). 1 lO'iaiiivrn ;ii i'''j iipi'asli •,h('ct occurs until the jets co deace At llii^ nouit, 
the anu.i-^li - hect (V i-^e^ to 0\iai Ikuico. the ncaanal or vertteal moinonium in the 
up ,1-^h i-^ !iu ion tl ^ anu .i mode' . • a^-t bo [)i ovuled 


The condition for no upwash flow, in the vertical plane of symmetry contain- 
ing the jet stagnation points, is presumed to Ije the nozzle height above ground wliere 
the jet boundaries touch at the ground effect height It is assumed that the 

U 

momentum per unit angle in the upwash sheet has a parametric form: 
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In the limit as s/d—^oo . fj^j^ — ^1. It is also assumed, that until the jets begin lO 
coalesce, the momentum per unit angle in the ujiwash sheet is conserved, or equal 
to the jet momentum and for any quadrant of the upwash sheet, 
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A function that satisfies all of the alxive conditions can be derived as. 
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Thus, ulicn tlie imiangement radius of the jet cijuals the half spaci ig of 
tiie two jets, the abo\e cfiuatinn \ieids a zero mor'entum flux along t!m ip u “ 

: rt c -mline. 

1 iguro 9-19 s/ious the variation in given b\ thiij model as a luielion of 
angle and llio parameter a^^. 'J'bus. this model ii wicatc=- lh.it as the sp. icing gets 
small . ‘1 or the nozzle height gets huger the momentum in the upv .isli sl'cct is ro- 
dislri 'iteri awa\ uom the ccnteiqxiint ljetv\een tiie two jets. The al.c-e cxpicssion 
L in be iuithei inleni.'.ttil to \u Id t''e noim.il momentum in the '.ipwasli sheet as- 
simune t!iU the upwash sticMiiih nes renrnn .ipju’o.xin atel\ rao’ial. 
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The model then i3redicts that the normal momcntuir in the entire upwash 
sheet decreases from a value of 2M^ to 4 , 

tl __ tJ 


At nozzle heights where the jets begin to merge, or their boundaries begin 
to ovej'lap, another m<3del must be hyi^othesized. In tliis situation, it is assumed 
that there is a region lx.-tween tlic jets where no upwash flow exists. It is no longer 
assumed that the momentum in the entire upwash sheet equals the momentum In one 
ot the jets. Hence, for r^^ > S, a function of the following form is assumed; 
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Tlie %ali e ot 2 at = m/2 has ixjcn held fixed and the angle v>u denotes 

o 

the angle at which upwash flow exists. This angle is estimated from whatever little 
empirical evidence that is n\n!lablc as, 
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and IS a function of tnc degTce of o\eiiap of the jet boundaries at the ground effect 


height as reflected in the coalescence jiarameter a . In the limit, when 

>> the angle ^ /2. The normal, or vertical U])\\ ash momentum, 

ill uvncral is ° 
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Thus, the vertical momentum in the upwash sheet vanishes rapidly as the 
jet boundaries overlap or as the jets coalesce as demonstrated in flgui'e 9-19. This 
tehavior is a result of the dlvetting of momentum flux into streamlines that become 
more parrallel to the ground. In the limit as v'y — b» ”/2, the amount of momentum 
in the upwash streamline does not have any effect*on the vertical upwash momentum 
because this streamline lies parallel to the ground. 

Further improvement in these simple models entails further empirical study. 
Unfortunately, It is extremely difficult. If not im.posslble, to measure the upwash 
momentum. Hence, the best ciddcnce for such models comes by way of correlation 
between ground up>vash pressures and total pressure in the upwash sheet. 

9 . 3 . -1 Ge ncivai Approach to the Pi'odiction of 'Avo - Jet Ground Pressure Distribu- 
tions for Close - Jet Spacing. - Figure 9-20 shows a schematic of the ground pressure 
phenomena that occurs for close jet spacing where the jet impingement and upivash 
deflection zone interact. Tlie upwash deflection zone now includes a portion where 
the upwash pi'cssurcs arc equal to the jet impingement zone pressures. Along this 
)-)Oi-tion of t’no ujnvash boundary the jet impingement pressures do not fully recova-. 
Hence, tlie Impu^jing flow does not ix-come eptuely parallel to the ground before it 
Ixgins- to deflect lipward and stagnate in the ujm'ash zone. 

To estimate the belunior of such a flow, it is assumed that the jet impinge- 
ment or deflection region is pezturbed due to the closeness of the tivo jets. The 
ambient pressure lioundary is assumed to ha\e the form: 

- e(l + cos(^) (9—18) 

The upwash Ixjundary 1? assumed to keep the same form except for the 
Ijcrtnrliation paramoter <r, or 
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Heiomng to figuie 9-20, it can Ix seen that Uie pi'cssui'e at the mutual 
boundaiw Ixtween the upwash zone and imjaiRgemcnt zone is abo\e ambient until tlie 
intersection of tiic i.vo boundaries where Uio picssure returns to .mibiont . Tlic 
jet iiiiiimgement r.idius to the upwnsli boundary is given by. 
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The reUitionshln Ixitwccn the h\o an}.'les and ^ is such the.t, 
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For a given ( and c7, the intersection of the two boundanes enn be found 
troni tl)e solution of the follow ing equations; 
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this gi\es ilu’ U\o a.ngles ^ and \\hich comnlctcU define tlie inteisection of 

uo 

the b\o lamiiul. lines. 


0. 3.-1. 1 .let HcAectK'n Hegio n. - To dc\ clop the anahsis, it is assumed that the 
uifogna.l ot the pressure in tlie jet 'PipiiipeineiU /one is Mil! equal to tiie jet thrust, 
or 
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Tile jncssuie intour.il is di\uled into b'.o jeuts. 'i'he first integi’a! renre- 
seiits the pn"-su; e in tiie region of intei action between the two deflection /.ones uhcie 
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a minimum pressure boundary exists. 
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Substituting Into eq. 9-52, 


V , (<P) 

: .. 0 .. min 






7T 1 > 

+ / /-> '■ )T? dl? ll(/5 , 

•'(P-’ Op gw W w ^ 


(9-53) 


IX't 


G {(fi) [ 

mi.'j Jn 


V . ((£>) 
nun 


r (</>) 

f (n )j) dT) where V {(p) 

g w w w nun r {<P) 

P 


1 


and 1 ^ f (I? ) »? dTj = C (n ) 
g g w w w S)lg g 


'1 lien the above expression reduces to, ujxin non dimenslonalization. 
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Unfortunately, the first part of this integral does not lend itself to analy- 
tical integration. Aside from its complexity, the implicit nature of - ' '’•‘tli re- 
spect to v^u makes numerical integration a necessity. 

The second integral in eq. 9-54 can be carried out and is equal to: 
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9. 3. 4. 2 Upu ash Do ricU ion Zone. - In the upnash deflection zone, the normal com- 
ixmont or \citical momentum in the uinvash sheet is assumed given by the momentum 
models liescrtbed m section 9.3..). 

Ilcncc, for a quadrani of the ui»n.'sh zone, 
r/-> \ 
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This integial is also diMikal into h\o iiai1;s, being 
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In the inteiaetio'i region of the uiwash /.one, the upwash pressures are .IS- 
C' limed gi\ ei! in , 
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Substituting Into eq. (9.5G), 
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Now, lotting 

Kq, (9-59) becomes. 


~K . ~ 0 C 

/. u 


AP 

uniL ’ , V 

=- ~ AP ,.(<;? ) 

AP. wall u 


-'ll 



r/v 


\ 


'0 cos''</> 


;n / 


\ 


/ 


AP (c7 ) -I f (n .1 C 

wall ^u \ / g ^ 


N , 


f (n ) 


d(P 


"/•J 


' i 




AP ,,('2’ Iciy ' 

.i .vail u u 1 


w lie re 


M .. ( Mj 


'' I 2 r 


(9-GO) 


*'-20 


Unfoitunatelj', these Integi'ais also do not lend thcmseh’es to analytical 
integration and must be numerically evaluated, 

9.3.4. 3 CKerall Solution. - To obtain the solution for the ground pressure patten for 
a given spacing and no 2 ./le height, the following two simultaneous equations are soUed 
iteratively for € under, 
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where and Ij., aie the two integiMl expressions for tlie jet inipingenicnt 
region and 1^^^ and coviespond to the integrals for the ujnvash '/.one. 

The Newton Itovative sclicmc co i&ists jf. 
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whcie a I'clavation t. it ten i j, ni.-.\ be nccessaiv to sneofi u)i con\ e’eciu'C. 

rimncs 0-il and 0-J'J .sliow ,i conqiansoii w itn test data of the computed 
giound I'l'cssiiic di.-'tnbiituin .tbine iho line id Ihc jcl st.'.geauon ptui.ts t.\ = i>) for 
\ n lous )Ct spatiMgs uui a u.'v '!e licight of '> d^. Heasoiiable goot! agreement is 
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achieved down to a spacing of hvo diameters. The assumed form of the jet 
deflection zone pressure function seems to differ somewhat from the test data 
but the maximum upwash pressures and widths seem, to be in good agreement. 

In addition, the predicted minimum pressures bet^veen the jet and upwash 
zone are in good agi*eement. Figure 9-23 shows some extreme situations for a 
close jet spacing of 2d^ and nozzle heights up to lOdj^. At approximately Sd^ 
the upwash momentum model predicts no upwash between the jets. As the nozzle 
height is increased, the upwash width decreases significantly at this spacing. The 
test data is in fair agreement with the predicted pressures. 

In general, the exponent of the upwash pround pressure function is 
assumed to lie 
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Further, the exponent assuma! to l>e a function of the minimum 


pressure in the deflection zone for 


To determine at forPmin>0, 
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an cniiiAalcnt profile for the upwasli pressure function is assim'.ed to exist whose 
integral force per unit ’ength is equal to the partial intcgiMl or partial lorce of the 
fulh de\ eloped pressure function (ougl'D = 1.5) 
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Giveii any I’min value along the apwash, equation 9-G4 can be solved to find 

the e'cponent a . When Pmin = 0 at = (i> , or at the intersection of the impinge- 

ug u uo 

ment and upwash deflection zone, equation (9-G4) reduces to an identity and, hence, 

cv = Q 
ug ugFD 

It can lx; seen tetli in the test data and in the predicted pressures that 
negligible interaction behveen the deflection zones occures at spacings greater 
than G dN. Significant interaction does not take place until a spacing of 3dN for this 
nozzle height as Is reflected by the minimum upvvash pressure level. 

Figures 9-24 and 9-25 show a sami>le of the entire ground isobar pattern 
for a lange of spacings and nozzle heights. Figui'c 9-2(J shows the kind of corelation 
achieved for one case. 

9.3.5 I'pwash Doca^ Model. - Figure 9-27 shows a schematic of the upwash sheet 
model and tlie characteristic scaling parameters. The upwash sheet is assiuned to 
be n reflection of the wall jot flow into the vertical plane of symmetry lying between 
the icts. Tile radial streamline jiattern of the wall layer is assumed to continue into 
the up.wash sheet and to a first approximaticii unpcrturlx-d by the turning region. 

'Ihe w.ill let tlow is assumed to enter the u)>w ash dcHection region with a cluirac— 
tciistie piofilc, half \elocitv wicHli am! maximum \ -'ocit\ or total pressure. The 
chai.u to istic length scale in the upv.a.sli sheet is the half \elocity width of the 
incident wall let ia_\er estimated at the wall location. Tlic piessure rc'co\ery 
region in the ujiv ash shoot is assumed to l )0 approxinritelv 3 times the half 
velocity width of the incident wall jet profile. The magnitude of the momentum 
flux in the ■■ esulting iijawasli stroaniline is assumed equal to that of the ineidcmt wall 
let. Ilcnce, the upwash sheet is ticated in a similar fashion as the wall jet with 
a tew exceptions. 

Due lu the high turbuloi.ee levels typically nioasurcd in the upwash sheet 
the fiow IS tonsideicd to I'o fiilh tuibulent and similar. Tne upwash \elocilv profile 
IS taken to l>e tl. it of a tree slu'-ii jirofiU" with no internal momentum defect due to 
tile a all hue.-, 'the upwash slu’ot giowth rate is assumed to be constant without 
ai)\ iiatia! tiaiisition phase. It .s .dso ao.-umml tl.at the m.igintudc ol the momentum 
fiiL< dmtiihution in the upw.i"-li slieel is gi\en b\ the iuiw iousl\ discussed momentum 
models. \r. c'xnrrssi.^n .malo;ous to that denied tor tne Wall let can now 'oe denied 
for the upw ,isli sheet: 
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R = .441 for a =1.5 
u u 


The velocity squared profile is assumed to be idenL cal in shape to the 
static piessure profile and as a result only one siiriilar'ty constant Is required: 
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The procedure for starting the upwash flow is similar to that of the wall jet 
whore an initial upwasli half width (bu,j^^) is determined from eq . 9-64 given an 
iiniscid starting maximum velocity. A motlel for the iipivash sheet growth must then 
be provided based upon empirical obseiwation. 
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The value of = . 30 is consistent with empirical obscrx'ations that the 
upnnsh flow dce.a\-s at a r.ate approximately tlircc times that of a free jet. 

An approximate expression c:in lx? derived for the upwash decay. Assume 

that the legion of intcicst is beyond ihe tiu’nmg region of the upwash sheet such that 

f = 0. Then cq. 'J-()o reduces to. 
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Norulimensionalizing boUi of those expressions by the value at 

<P = 0 or X =0: 
u u 



The more general procedui'e was used to generate comoax’isons with test data. 
Figxirc 9-28 shows the upwash ccntcrlrne (X^ - 0) decay and half velocity boundary 
growUi characteristics Ln comparison witli some test data. The model predicts higher 
initial total pressures for tlie closer jet spacings along with a smaller tuimhig region. 
The decay curves tcixl to merge witli increasing height above ground. The model half 
velocity boundary growth stalls off at a smaller value closer to the ground for closer 
spacings and then grows at the same ra*e. 

Figure 9-29 shows the predicted lateral variation in the upwash pressures for 
H/dn = 4.0, S/dn = G. 0 for several heights in the upwash sheet. The total pressure in 
♦ lie uiiwash siicct drops off more rapidly compared to the centerline value for greater 
heights due to tlic radial spreading of the upwash streamlines. Figui'c 9-30 shows r. 
compnnsoii of the Litoral total pressure decay with some test data at h/d>g = 3, 
S/dvj = 4. Figure 9-31 demonstrates the small influence that nozzle height above- 
ground has on the laieral cliaraoleristics of the upwash sheet ai a fixed o’.ane. 
1 igure 9-32 shows the lateral decay characteristics for a close jet spacing of 
2dy just jinor to significant coalescence effects. At greater hciglits above 
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ground significant upwrah flow does not exist along the centerline between the 
}ets. Figure 9-33 denicnstrutes this effect where the up^vash total pressure now 
has to be nondimensionalized by the jet stagnation point pressure. Figure 9-34 
shows an even more extreme condition with some test data plotted for comparison. 

9.4 Model for the Prediction of Upwash Forces on Slender Bodies 

it Is a well known fact that when an upwash flow field impacts on the underside 
of an aircraft an unw-ird lift force can be generat'^d on the aircraft. This upwash lut 
force is soinetinits negated at heights close to the ground by negative interference 
forces (suckdown) caused by the secondary induced flow field gcneiatcd by the 
entrainment action of the primai’y jet and wall jet flows. It is the intent of this section 
to dovc'op a method Cor the prediction of the up.vash lift on a slender fuselage where 
the suck(lou7) force has a minor effect in compai'ison to the upwash force. 

9.4, 1 Upper limit of Upwash Forces irom a Momentum Capture Model. Figure 9-30 
shows a «'kctch of the flow model to be used for momentum capture; 

Tlic lic'iglit above gi'ound of the fusehige and nozzles aic assumed equal. It 
will be assumed that the entire vertical momentum included in the upwash sncct by the 
length of the fuselage can be converted into a lift force. Accordingly, this lift force 
car/ be comouted as: 




(a + 2(1 -a )sin <p ] cos cp do? 
u u u u u 



(9-72) 
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In tlic above expression, the lift force nas been nondimens ionalized by the thxaist of the 
two jets. Hence, the maximum lid force attainable for an infinite length body with 
au = 1 is 



This expression, in the limit as^ — ► oe tells us that only 31. of the 

X i 

thrust can be achieved as the optimum ideal liftforce neglecting jeteoalescence effects. In 
reality, for a finite length body, the lift for*ce will be markedly less than this optimum 
value. Eq. (9-72) is plotted versus nozzle height for several jet spacings in figure 
9-37 for a t^iJical T>pe A fan jet V/STOL fuselage length of 8 d^.. The jets are also 
placed at the midpoint of the fuselage. The figure shows that the maximum uowash 
force on this finite length lx>dy is about .25 Tj, Tlie upwash force also dropi off 
rapidly for a closely spaced jet corJiguration due to the rapid coalescence of the jets. 
The momentum captiu’c model neglects any finite width or fuselage shape effects. 

This t\pe of calculation is most representative of a fuselage fitted with 
strakes which is capable of approaching the force as predicted by a momentum capture 
model. 

1) . 4 . 2 General Appi-oach to the rrecliction of Upwash rorcc on Arbitrary Slender 
B odies, In general, the fmiio w'idth of a fuselage will have a marked effect on the 
amount of momentum iransfcrred to the body as a lift forae. In addition, the shape 
of tltc fuselage cross section will determine the local pressure distrib'ition around 
tlic bcxl}’ and thus also effect tlie drag or upwash force. A general approach to a 
prediction model for tlie upwash force on slender bodies is developed m this section. 

9.4.2. 1 Maximum Stagnat ion Pressure D istmbutlon Along Centerline of Fuselage 
Underbodv. Figure 9-38 shows a schci.iatic cf the fuselage, the pertinent coordinate 
svstems. and geometrical parameters. The centerpoint or centerline of the body is 
locatetl at an arintrarv liciglit Zp reiatne to the ground and the nozzles arc located 
at a height h above the ground. The fuselage depth is given by the equation 

Zp-f(Xj,) 

The fiLSclage hrdf width is also defined as 

3 ^ = S (Xp) 
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The nozzles are placed at the midpoint of the body. The angle in the ui^wash 
sheet in:ludcd by the half length of the fuselage is then given by 


= tan 


-1 


£ 




The unf'orslde fuselage centerline coordinate relative to the upwash sheet at 
any p sitiori is given by 


Z 

u 


- S. 


+ Z -Z (X ) 
p B B 


To determine the maximum pressure distribution on the underside of ihe 
body, it IS assumed that the local stagnation pressure on the body is representative 
of tlie stagnation of the cross flow component of the ufAvash velocity. The cross iloiv 
upwash velocity is just the ncmial component of the upwash relative to the undorsido 
of tlie body. 
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y. -1.2,2 Cicucral Force integral. If F\i (X) l. defined to be the upsvash force per unit 
Iciigth ef ihe bod\ , tlicn (lie entire force on the l>-j'ly can be defined as 

f , 
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F. 


T 


f, ■ 


(3-7 1) 


0-3S 




\ 


Converting this integral to the upwash coordinate eystem entails, 

X„ 


Thus, 


tan (p 
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Hence, the integral eq. (9-74) for the total force becomes. 






F (o’ ) sec**o do’ 
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Now, it Ls assumed diat the force on any cross section is proportional to the local 
cross flow stagnatio J pi\!ssurc, or 


F (<fi ) = I’AP (O' 
u uB uo 


11 B^ j 


B./^ub’ 


G (\ ) dv 
P' 
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wiicrc Gp (y) is the pressure function rcpresentbig the local flow around the cross 
section. 


Uvxin substitution of this expression into cq. (9--75) and using cq, (9-73) for 
the stagnaiicn pressure, cq, (9-75) becomes. 
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If the integral over the body width Is Identified as the drag, and a 
sectional drag coefficient is hitroduccd as. 
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and cq. (9-77) is nondimensionalizcd by the thrust of the tv;o jets, the final expression 
for tlic total force Ixcomcs , 
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where die above lengtlis arc noiKlimcnsional. 

For a cyiindvical bcxly (i.o. , constant width and (iepd)) wiiosc imdcrcido is located at 
the nosczle height, or 


Z 

P 



= h 


tlio above expression reduces to, 
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9. t.2.3 Estimatinn of Liie Drag For Uiffcix'nt Cross Sectional Slinpcs. Figure 9-39 
shows appi'ox.iniatc drag models corresponding to a circular and rectangular shaped 
cross sci-tion body. The nioticis for the drag for the upwasii proolcni arc markedly 
different than conventional aerodynamics. A c>linuor immersed in an infuutc sticam 
will exhibit 2010 drag iX tJio flow is attached ami mvisoid. The \uscous flow. 
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on the other hand, will separate at some point aixl a base pressure will result causing 
a drag on tiie cylinder, hi the upwash problem, tlic viscous flow will impact the body, 
flow aroimd it and also separate at some point. The difference being that the pressures 
aft of the separate point will be the result of secondary induced flow. Hence, resulting 
in base pressures much closer to ambient as compared ^^o the body immcrsecl in an 
already energized infinite stream. These considerations were the basis for the drag 
models. This is also the primai’y reason for consideration of slender bodies because 
these small secondary Induced pressures acting over a large area will result in a 
significant suckdown force. For the present time they arc neglected. 

Figure 9-39 b shows the model for a flat bottom body. In this model, a 
finite width upwash flow impinges on the umlcrsidc of the fuselage. The model assumes 
that the flow separates at tlie bottom comer of the fuselage and continues without 
ciuwaturc away from the b<xly. In reality, the flow would have a tendency to follow 
the body after separation. The force can then be easily estimated as the portion of the 
force tliat would be exerted on an infinite flat plate but taking into account tho finite 
width of the underbody. 

If tho pressure recovery zone for the stagnating flow is estimated to bo 
approximately, 

y - a bull where a ~ 3. G 
o o o 

in terms of the half velocity width of the incident upwash stream. 

r o ,1 

Now, let C. (n)-'ll-n I 

‘ FLAT 

wliorc n = and u ~ 1.5. 

‘ 0 '^o^ull 

Substituting these oxprcs.sions into cq. (9-77) jiclds for the ci'oss-secticn.al 
di-ag fiuietion, 
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where - 


^ 1- 


The drag coefficient for the flat bottom fuselage is thus totally dependent 
upon the relative size of the Incident upwash flow arxl body width. In the limit of 
*^Uj^ — ► 09 or — *-0, the maximum drag coefficient is 1. Tills can be interpreted 
as bi the limit the stagnation prescurc acts over the entire width of the body 
(analogous to a pitot tube). ,0n the other hand, the minimum drag coefficent will bo 
obtained when Vp ~ 1 and yields a value of , 3 15f>; The drag coefficient behavior 
versus the piiramcler i)g is shown in figure 9-40. 


For the cirular ci’oss section, the potentLal flow 'solution for the prcssimc 
dlstriliution ai'ound the cylinticr Immersed in an Infinite st?'cam is use<l. The flow 
is assumed to separate at some body angle Gg and the base pressures arc taken 
to be ambient. Thus, for the circular cylinder. 
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It is .assumed that tlic flow separates aniwhcrc between IzO" and IGO". 

The drag coefficient as a fimction of separation angle is shoum in figure 9-41. The 

niaximum value of Die drag coefficient for a cii'cular cross section occurs when the 

flow separates at G = 150° and is C_ = .333. 

S DrC - - 

To see how well these model., work, some correlation with test dtUn was 

carried out for c\lindrical bodies w’hosc underbotlv was places! at the no:.zle iicight, 

Fij-aue 9-42 shows the correlation achieved for two such bodies for c ircular and 

rectangular shapes and a cross sectional width ecpial toono diameter or half width of 

one noz/le ladius. I'or tlic simple nioc'cls used, the correlation acliencd is quite 

reniarkablo especially smee the sue kdown' effect was ncglecteil. The test data drops 
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off at a slightly higher rate tlian that prallctcd probably because of the sucktiowii 
effect. Thus, tbe mysterious sensitivity to body cui'vativc obsei"V'ed in past 
experiments can now bo logically explained. 

Figure 9-43 shows some parametric curves generated for tlio flat 
bottom cylindrical body as a fimction of botiy width. As expected, the predicted 
force increases propoitionatcly to the width. Unfortunately, this treixl is an Ideal 
one .oiul at some body width tlio suckdown forces will predominate. Figure 
9-44 further shows tlic predicted trend for a parabolic body of maximum cross 
section equal to th.at of a cylindrical flat bottom Ixxly. As expected, tlio foice on the 
paralJoUc body is somewhat less tlian that on the cylindrical boily. 

These mo<tols could further bo improvcxl by Includmg tlio suckdown effect 
mid also extending the estimate of the drag coefficient to arbitrary cross sectional 
sliapob. 










I- igiue 9-4, - Flow behavior just prior to coalescence. 
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10. CONCLUSIONS 

In gi’oimcl effect twin jet tests were conducted to determine the gc.'eral 
features of the resulting flow field and to ascertain whether such flows arc amenable 
to a mathematical model. Some of the general features looked for were instabilities, 
basic flow structvurc as the jet spacing was decreased or height increased (causing 
jet merging), t^TiOs of moasuremonts , the deiioction zone which results from the 
collision of wall jets and the rcsultmg up'A'ash sheet, and the eff^'cts of bodies placed 
in tlic upwash sheet. 

The isolated t\vin jot tests r. vcalcd unstable up\vash flow over a narrow 
range of jet spacing in the vicinity of three jet diameters and jet heights around i 1/2 
to 3 diameters, bistabilities were observed in pressure measurements on the ground 
plane in upwash pressures, and in force measuiemeius on bodies placed in this flow 
field. 

Regions of below ambient pressure wore noted in many areas of the isolated 
two ict flow field. These included Ixith static pressure measurements made by tra- 
versing proiics across the upwash sheet and measuring surface pressures from the 
stagnniion line on tlie ground plane, fhe low pressures in the upwash sheet apocar to 
I’l a clmracteristic of 'he two joi ground imi>i.Tgement flow and were obtained for all 
valuc.s of height and spacing that were tested (from 2-6 diameters). The lowest static 
prc.ssures occur in the center of die upwash sheet and not at the edges. 

Three tvpes of borlies placed in the resulting two jet upwash sheets as well as a 
l/)0tli scale aiiviaft iiiodcl witli strakes of vtirious depth produced mai-kedly different 
)ct induced forces. Net litt forces over 20'. of the jet thrust were measui'cd due to the 
upwash sheet ui'puiging on fuselagc-t\ pc surfaces, Stiakes were bcnericial on simulated 
fuselage surfaces hut not on lai go flat plate surfaces. 

A tl- • 'ough im ostigcUion wa.s made of the flow field that results wlicn the 
ids arc cKisoh spaced and mciguvt oJ tlie outer boiuidaries occurs prior tc unpingc- 
iiien! on (he ground plane. It w.a.s’ concluded (hat the upv/ash .sheet continues to be 
formed even ai close spacings and dial i! sepaTaie peaks appear in the free jet , pressure 
profiles vviiho'al a groi.nd plane, so,'«.u aic impingement and ujiwash fo-mation will occur 
with a ground plane. 
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A new ficrtv model was developed for close jet spacings which accounts for 
non fully formed wall jets prior to the development of an upp.vash sheet. It includes a 
jet impingement zone wall-jet transition region which adequately predicts the effects 
of jet spacing and height above groiaid on the flow field pressures and upwash flow 
field. Integration of the upw'ash sheet pressures were used to compare with the test 
data on fountain lift. Good correlation was obtained with simple twc-dLmcr.slonal 
bodies of varying cross section. 
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Appendix A 

Computer Program Description 

Program 

Name : GRUMJET2 

Purpose: VTOL Two - Jet Impingement Interaction Probiem 

The program is designed 1o estimate the flow characteristics associated 
with two vertically impinging and equal strength incompressible jets. This pro- 
gram is specifically oriented towards the closely spaced jet interaction problem 
where the deflection regions interact until eventually jet coalescence occurs. 

Aside from the basic flow characteristics, the program assumes the 
symmetrical placement of a slender fuselage in the upwash flow. The upwash 
lift force is then computed for u cylmdi’ical fuselage of constant cross-sectional 
shape. The body parameters, in terms of width and depth, do not vary longitu- 
dinally. The program estimates the force baseci on a rectangular and circular 
cross section. Two planes of symmetry ore assumed and sill output applies to 
one quadrant of the flow field (ie. equal jots and nozzles located sit the midpoint 
of the fuselage). Jet entrainment effects may be significant but are neglected 
in this program. Hence, the force is only that due to upwash impingement. Some 
residual prog'rnmming exists in the c*ode for a parabolic body of revolution. 

These cards have been commented out but may be used if desired. 

Input Description 

Note: All input parameters are in terms of nozzle diameters. 

Card No. Code Names Format 

1 HD. SD, ZPLD, DZPL, ZFINAL 5ri0.5 

Name Definit ion 

HD Nozzle height above ground 

SD Nozzle spacing 

« 

Note The program will compute one or sevcrnl positions of the fuselage relative 
to tlic ground foi a fixeu nozzle neight nnovc ground. 
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ZPLD 


Initial fuselage height above ground 
DZPL Increment in fuselage AZ above ground 

ZFINAL Final Z coordinate of fuselage height relative to ground plane. 


Card No. 
2 

Names 

XL2 

WCON 

ZCON 


Code Names 
XL2, WCON. ZCON 
Definition 
Fuselage Length 
Fuselage Width 
Fuselage Depth 


Format 

3F10.5 


' Note : ZCON determines position of fuselage underside rel.'’tivc to it: ZPLD 

locution. Bottom of fuselage will be located at ZPLD-ZCON at first com- 
puted location. The up wash sheet properties are also computed at this 
Z location. 


Card No. 


Code Names 


Format 


3 


IPBAR 


II 


Name Definition 

IPBAK Integer controlling the output of ground pressure 

pattern. IPBAP - 0, no jircssure pattern output 
= 1, pre-^suro pattern output is desired 


NoUi_ Card No, 4 is not required if IFBAH = 0. 

Card No. Code Name Format 


1 

Name 


NU 


Nu 12 

Definition 

Number of pressure values to be input for computa- 
tion of ground isobar pattern NU < 25 


Note Curd No. 5 is repeated NU limes. 
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Cord No. 


Coce Name 


Format 


5 PU F10.5 

Name Definition 

PU Nondimensional pressure lor isobar pattern 

0.0< PU< 1.0 

Figure .\-l shows a typical input set. 

Printed Output Description 

Figure A-2 shows a typical computer printout. Most of the geometrical 
output quantities are nont'.mensionali/ed by the nr^zle radius. All velocities 
and pressure arc initially nondimensionalized by the nozzle exit velocity and 
stagnation pressure. All pressures aic relative to ambient conditions. 

Note: RN and D refer to nozzle exit radius and diameter. VN refers to nozzle 
exit velocity. 

Input Parameters 

The first set of output echoes the input parameters. 

Output Titles Defi’^itions 


11/D 

Nozzle height abov'e ground 

S/D 

Nozzle spacing uistance between jet centerlines 

Z/D 

Initial fu.=elugc height 

DZ/D 

Increment in fuselage position 

ZFINAL/D 

Final location of fuselage 

L/D 

Body length 

W/D 

Body width 

ZB/D 

Location of underside of body relative to specified fuse- 
lage location 


Note 


If /.B - 0, underside location is coincident with fpecified fuselage pnsiticn. 


IVVTtC< 


Jot Decay Rofrlon . 
Output Titles 

di;lg/d 

ZPC/RN 

ZFD/RN 

Z/RN 

RJH/RN 

RJ/RN 

ALP 

CV2 

VJ/VN 

Ut:OHE/RN 


Definitions 

Jet ground effect height relative to ground plane 
Length of potential core 

Length of potential core and transition i*egions 

Jet axial location measured from nozzle ox't 

Half-velocity ra«lius of jet measured from jet center- 
line 

Half width of jet 
Exponent of velocity profile 
PZomentum flux intogial 
Centerline velocity 
Potential core radius 


U't Oetlcction Region 


Single ]et iii.pingcnient chaructcnstios arc printed. 


Output Title s 
UOll/UN 
RG/KN 
RO/RN 
VG/VN 
DPS/DPTJ 

•\r PG 
H 'UGH 


DPS /DR I'G 


D efinitions 

Half velocity radius of jot at gi'ound effect height 

Jet half width at ground effect height 

Dcfiection zone or pressure recovery radius 

Square root of ict i''round stagnation pressure 

Jet ground stagnation pressure nondiniensionahzed by 
nozzle stugnotiou piessure 

Exponent of ground pressure disirioution 

Radial location in lienection region leiative to ground 
stagnation point. Komhincnsionalizcd by groimd eflcct 
half velocity radius of jet 

Ground pressure lu'ndimcns’unahzeci by gtounci stag- 
nation pres.sui ‘0 




Output Titles 


DcfiiiitioMs 


R/HN 

Dl’S/DPTJ 

Wall Jet Rojrion 
Isolatud 
Out ! uit TUlos 

i)i:ls/un 

VM/VO 

AU’WO 
JtWOH/UN 
inVO/'KN 
R /ItN 

VM/VC 

VM/VN 

mvil HN 
H\V 'HN 
Dn.iU./ltN 
1)1 1 r.) 

KDl I. 

VI'O'M it 


Hiidiii) locatk»n in defection re^pon uor.dimonsiouidiiod 
by noi/.lo radius 

Ground prossura nondiincnsioiiali/.cd l>y iiO//.lc staj'ua- 
tion prossuro 

wail jot prupoi'tics arc printed in th./. set of outp-ut. 

De finitions 

IJoundary layer tluekiies.s at stapimliou point 

Inviscid iiia.v.iiium velo*,ity at jet luilf uidtli (lUill) 
radial location or start of tran.^ition to turbulent 
wall jet 

Exponent of velocity profile at start of wait jet 
Hall vciocilv thickness at start of v;a)l jet 
Initial thickness of wall layer 

Radial location o'l wall ;el measured fiMai jtroinul 
staitnation jiuint 

Mavnauni veUicity in wall iaver nondira’iisioiiali /cd bv 
square rot^>t of pronnd stu^jiiatuai ju'essure 

Ma.Miiuim velocity in wall layer referenced to iio//,le 
exit xe'ocity 

Half \elocity thickness of wall laver 

I'liickiu’ss ot wall laver 

Houmlarv la\a'r thukness 

Maxi'iiuni siapn.slion pressure in wall laver 

iiatio oi boundary layr ti> total \.a!l layer llneknes-, 

1 \poiieii! of bound. irv lav> r ju'olile 


ll-l. 


Two Jet Impiiipfement Interaction Output 


Maximum Ground Pressures olorifr Upwash Ground Stagnation I.ine 


Outpux Ti tles 

XW/:tN| 

XW/S I 

PMAX/PJ 

PMAX/PMAXO 


Definitions 

Coordinate along stagnation line measured from up- 
wash stagnation point 

Stagnation line pressure nondimcnsionalized by nozzle 
stagnation pressure 

Stagnation line pressure nondimcnsionalizod by up- 
wash stagnation point pressure 


Upwash Momentum Function 


Output Titles 


Definitions 


HG/ilN 


Radius of jet at ground effect height 


AGON 

XMOMZ 

Pino 


Constant in momentum function 

Total vertical momentum in upwash sheet nonuimen- 
sionahzcd by the optimum value of Mj/2n. 

Coalescence angle used in upwash momentum model 
If PH 10 = 0 jets are not coalesced 
!f riilO > 0 jets have begun 'o coalesce 


Note : Depending on the value of PHIO (i.e. zero or non-zero) the constant RCON 
appbes to the appropriate upwash momentum model. 


Cll Nondimensional upwash deflection region width constant 


Note : There are two possible outputs that can occur at this point. If the jets 

• ire .spaced far enougli apart the comment: 

.IL'T A.N’.':) UPV.ASH n.'.rLFCTIO.V REGIO.NS DO NOT 
JN'T/ R.iCT. CL’ F.STJMATH IS CORRECT 


111 t'ua case ti5C •. .'-I'lc c! ClI .s correct and the ict laipingcment and up- 
deflection regions lU'C independem. The perturbation parameters 
.irc then definid a-, ; 
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EPS = 0.0 


SIG = 1.0 
PHIUO = 0.0 
PHIO = 0.0 


PMIN = 0.0 alonp: upwn&h line 
ALPUG = 1.50 

If the nozzle spacing and height above ground arc such that the deflec- 
tion regions interact, the subroutine INTERG will be called and 

♦♦♦CALL INTERG ♦♦" 


will ho displayed. The following output will be printed. 

Output Titles DMinitions 

Iteration cycle Number of iterations required to find solution using 

Newton's method 


EPS or El’I 


Perturbation parameter for jet impiniremcnt pressure- 
distribution 


SIGMA or SIG 


Perturbation parameter for upwash doncction region 
presburo distribution 


CU or CL’I 


Upv/ash width cstiraato prior to iterative bolution 


PMIIJO Anglcb defining the intorsoction of jet and upwash 

PHIO deflection repons 

Noio The following values apply along the line c«jnnccting the jet stagnation 
points on the ground (i.e. x = 0). 


DHLPWO 


Pressure at upwash stagnation jjoint nondimcnsionahzed 
by nozzle stagnation pressure 


PMIN Minimum pressure between jot and upv ash aoflcction 

regions noridimcnsionahzod bv na//h: stagnatiuii 
pi‘Cbbure 


AI.l'LUJ 


Lxi/ijuent of upwash ground pressure distribution 
function 


(CL’‘) X (SIGMA) I-'inal value of upwash thickness eonstam 
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Note : If IPMAR 4 0, the followinfj output will occur. 

ComputdtioM of Two - Jet Ground Iso»u»r Piittern 


Outuut Titles 


PBAR 


Definitions 


Input values of pressure, nonriimensionaliicd by 
nozzle stuj'nation pressures, will be echoed. 

Note : The followint? will be repeated NU times. 

Output Titles Definitions 

IJKT .Number of points on jet iniping^onient rej^on isobar 

ID Number of points on upwash deflection rcfpon isobar 

If IJllT - 0 or ill := 0, the .spocifu'd value of PDAR was not found in ground 
pressure distribution 

Note ; The foltowinp ocxu’dinatos arc referenced to the jet .^tairnation point on 
the (Tround. 

XISOJ X co<inlinatc of jet isobar 

1 IsOJ V (.'oordinate of jet isobar 

\lSOti X coordinate upwash reppon isobar 

Y1M)1’ Y coordmalo of ujiv.’ash repnon I'-obar 

Ihe final set of ouunit in this section is the upwash deflection 2one line. If the 
dellcction zones do not interact, aaibicnt conditions exist ulonji this line. 

M'l’ X coordinate of upwash line 

Vl'P Y coordinate of upwash line 

l.’cnouiation of Gjiwash 1 low I'lold 

I'p'.sash Streainhno Projicrties 

Note l '.\o streamlines are printed 

I'll! (Mil Titles Defimtioiis 


IMIli) 


!\\ . \ 1 1 ,. 


A/iiiuithal anple of u’pvash streemli'ie relere'iced to 
jet ijrouiu! c»K.>rdinato system 

l:i\'isc;d Mirr.nu' .’Otrion maxinium \clocitv 


1 1 -;i 


DELPVVJ 


BWGH/RN 

ZOU/RN 

BUO/RN 

BUOHN/RN 


Maximum pressure on upv/ash stagnation line on the 
ground where upwash streamline orginated 

Half velocity width of incident wall jet streamline 

Upwash turning region height above g’ound 

Initial upwash width 

Initial upwash half velocity width 


Upwash Streamline Decay Properties 


2/RN 

VMU/VN 

lUJM/RN 

BU/KN 

DHLPU/DELPJ 


Upwash streamline coordinate (ZS) 

Upwash maximunj or centerline velocity 

Half velocity width 

Half-width 

Upwash centerline or maximum total pressure nondi- 
mensionolized by nozzle stagnation pressure 


Upwash Properties Computed At Z Location of Underside of Body 

If ZB = 0, the output will yield the upwash properties at Z = constant 
plane above ground. 


Out put Titles 
ZU/RN 
RW/HN 

PMIB 


X/RN’ 


Definitions 

coordinate measured fixmi ground plane 

Radial coordinate from jet ground stagnation point to 
upwash stagnation line 

Azimuthal location of upwash streamline referenced to 
jet coordinate system 

coordinate in upwash sheet or X coordinate on 
fuselage 


VMU/VN 


Maximum or centerline upwash velocity 


PU/PN 


BLi.’/KN 


Maximum or centerline upwash total pressure noi'di- 
monsionalized by nozzle stagnation pressure 

Half- velocity width of upwash 
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Upwash total pressura nondimensionalized by total 
pressure on streair-Une oripinating from upwash stap- 
nation point (X = 0, Y = 0) 

Upwash total pressure nondimensionalized by jet 
pround stapnation pressure 

Stapnation pressure on underside of fuselape placed 
in upwash sheet 

Upwash lift Force 

Two values are printed: 

FLAT BOTTOM VEHICLE WITH SHARP CORNERS 

BODY WITH CIRCULAR CROSS SECTION 

Output Title Dcfimtioii 

I.U/2T.T Upwash lift force nondimonsionnlized by the total 

thrust of the two jets 

Computer Profn'am Listinp 

Fipurc A- 3 shows a fortran listinp of the computer program. 


PU/PUO 


PU/PJS 

PBS/PN 
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Fifpare A-1, Sample input Data Set 
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I'OO I UN J I NIM 

X 1 tUf. - ( 111 ••*HMIM <' 00 *r 
‘UIM I = O . 

Nil! TSM-NUPT‘ I 
on Sv>0 I I 5 I • N IK T ',M 

KOI f Kll *1 f I <- < 1 . K (l!K r H I < 1 U4 1 > 1 I 
»?{)f f f «0 * 

M1» M - MM f I : 0 t 1 . K'UM f M I I I I ) ) 7 
KOI MSU Mil MOO J 
f M (vUI M'.UOi ‘.HI ‘ > 

If Kill f J > I MI' 

Of MI I HI I 1 M 1 ' PMl I I J I 
Mini MiMlf GOOIMIOMMMI } 

<00 I UNI I NUf 

XII U. Ot L t GOMIM I 

r.M I S I M ( O . O * Al I 1. » K A I .Ki • ( L I lu ) 

X I ' •!. ( sr '00 01 I f ‘.0 < < f I » »U I > A ( Ml« I Ml« . Of • I > « 1 . \ «| I I »«.* ) 

I *1 I I • \ KIM M I ' . I NO MI'M ”.0‘, T Nl 7 . »» HUM *t » 1 **.' > 

X I I ' V 1 t Jli f X I lU 

f K < I I K ) ^ I 1 IT 

I»l M I U t . *,$» I I M 1 HO ) '\f 1 S 

Kiinu.’ o c 

I M f M I I .U 

V t ill I (. 1 « • 1 N ( * iM un ' > 

I M I ( 1 Mx ) '• 1 1 I'll M I 'U< . X : An 

1 » » M i At. I u t ' It' 1 U '» C> / 

* I I I K Mi t > l.' I Ml 1 OOP 

M' *1* 

OM 1 I A» M 

O- HI I .M n 


I i‘un u \ -.'i, ‘/r.nlnuir‘1 


1 NI Ol/A to 
INIOOA 10 
INI Ol'AGO 
IN 100 AA 0 
INTOOA/O 
ISTOOAOO 
INTOOAVO 
INI 00/00 
INIOO/IO 
INrO0/70 
INTOO/ to 
INt OU/ Av> 
INI 00 /GO 
I NfOO /AO 
INI OOw’ / O 
IN I OO^ .TO 
I N I Ov' /VO 
I H \ '■'OOOO 
IN) OOOIO 
IHTOOUJO 
I N I OOIMO 
I NT OOU40 
INT 00H*.»O 
INT OOUnO 
INT ooa 'C> 
I N T OOUGO 
I N I OOOVO 
INI OOVOO 
IN TOOV 10 
I N I v'v •> .‘O 
INT ► G> 
I 'O 0 ; 4 <> 
IN I 

I N T OU“ . ) 
IN TOOV V 
I N I OOVhO 
I N T v>0 /VO 
INTO! OOO 
I N . O! O 1 O 

1 N T <M ■>:*o 

1 N 1 O T O 

IN T O 1 040 
INT 0| OGO 
X N T O I OAO 


1 N I (' 1 

O ’O 

IN ! O 1 OHO 

I N f O 1 0'» O 

T ffl O 1 

1 no 

INTO! 

1 1 O 

1 N 1 •' 1 

I .’O 

1 N 1 <M 

1 ti# 

I N I <> 1 

1 10 

I rM 0 1 

1 • O 

1 N UM 

• .O 

1 N M' 1 

1 ’O 

1 N 1 *' 1 

1 U) 

i N 1 O ’ 

1 • *' 

I r M 1 

u' 

M« 1 > 1 

. V Kt 

MM v> 1 

o 

M' I O 1 

. * o 

M 1 •' 1 

4 0 

MM •' I 

GO 

1 N 1 O 1 

.v> 

1 N . . 

’ / 

I N M ' 1 

'K. 



1 


< XTR.CO. ?> rUIE-« <rui < I TK > - riJI < MR - I > >/DL*f I 
IF <ITR,EQ,4) rule =<rui ( I7K >-FUI ( ITR-1 J )/OCUI 
IF <iTR.co,i»> f rr-u Ji < 2 IK ) -iji < I iK-i > 

IF (IYR.Ca.4> r Ml <rjl < I TFx >- f ) l < IlK i>)'DEft 
2000 CONTINUE 

ttCT -ruii.wi Jic* r jiEe»Fuic 

<tji( i>«Fuic-rui< i >«rjic)/Pf:T 
MELCriFui < I )*f JIE r JX < 1 )*ruiL )/£il.T 
EPlJ4l>-'CP< l> ♦RFL 
c<j4i)«c( ») 4 ki:l AK*nnt-C 

VO*. rURHAT < lUX. ' I ICKAI UIN C YCt L . I i * 2 X • ' i_f S « ' • F I O • ts * 2X * * S I OMA« ' • ! 

ir < nri.r ) . LT . I .'t -3 . ANP . AfiS< fJELO .1. T . t .F -S> IFLAGal 
ir < If LA(J . CO. 1 » URIIC <A*90V) 

V09 FORMAT < JOX • ' ^OLUl 1 ON OF UROUNP PRC^^.GURE fi I S T R 1 HUT I ON HAS PECN 


IF < IF I AG. I O. I ) UR2 1L <6»90* 
900 CUN 11 NHL 
ST(JF* 

VO^ S1G»C<J) 

FF S=-Ff< J ) 
f‘Ull-«S I04CH1 

UM TF 1 1 ) cini 

fflKMAl < .'OX f ' U» UIA«'.H TMirKKTl 


J» rP( JF I > .C< M I > 


I r < I F 
I > 

KF 1 UKN 
f Nil 


MAI < .'OX f ' U» UIA«'.H TMIFKKTl.S LONATANT# <CU» X < S UlMA ) w ' , » 1 O . ^ 
< IFPAFx ,Nt .O ) LAI I GFI (U < Al f (. . PCI f S • SD . RO » I.PS • CU I »G1G»FHXU0» 


INTOirvO 
INTO J 300 
INTOl JtO 
INTO! \^'0 
IN10133U 
INTO I 5^0 
INTOl JUO 
1 NTO I 3A0 
IfMOl 170 
iNTOi ino 
4 10. INTO I 3VO 
INI 01400 
X NT OI 4 1 O 
INIOl 41'0 
F0U4NT01430 
I NT O 1 4 40 
INT 014:,0 
1 N1 A I 4<*>0 
X N T O ! 4 O 
1 NT(»1 4HO 
IN101470 
IHTOILOO 
I N 1 0 I 1 t> 
> 1N101G20 

r HIO I N I O 1 :,30 
I NT O 1 ^*40 
JNTOILGO 
IN! 0 t ^60 


^ Opt. I lied 





uuu 


StmRQUT INr G^'LOT < ALI^f* • DETU • SP # «0 # KPS t GUI .SXGl •FHlUO.PHXO> OPLOOOXO 

OtHCNSSON rU( pRU/HLL< 102) tOP) 0 Xf'L0T0(202>^ Y)U.07tJf PPL 00020 

1202 > tXPLOT J ( 202 > • YFLOT J< 202 > # XAX X S < 30 > » V AX 1 G < JO ) *XDATAC XO> » YDATA( ICPLOCOJO 
20 ) 0 PL 00040 

COhMON /MCXOHl/ MDflPPAR 0 PL 0 C 050 

DPLOOOAO 

THIS KOUTIHi: COHJMJTES THE UROUND ISOBAR PATTERN FOR THE TUO-JET GPL00O70 

IhPXKGrrXRT flow fiilp gplooomo 

cr L ooovo 

KtTA<CTA»A>»<l.'E I A$<A>4«4 OF LOO JOO 

GPLOOl (O 

INrUT NUMBER OK ISOBAR VALUl S CT LAO 120 

Of*LOCJ JO 

Rt02(fi»M) Nil GPL00140 

51 FOKMAl <12> 0»t00X!*O 

WRITE HU GPL 00 1 AO 

53 rORMAT < I ' . //JOX • ' COMPUTATinN OF TUO-JET P.ROUND XGOl*AR PAT T L GPI OO 1 70 


IRN '///J2X* !2i 

2// > 


VALUES or PRESSURE GrrClFIED FJR fATTfRN' 


lOBAK VALULS FOR CiRGUMfi PATTERN 


KIAli(5*l> (I IM 1 > • J « 1 f Nti) 

1 I CiKMAT I F t O , S ) 

W^^IL <M/< X > # I *1 f VU> 

112 inUMAI < lOXf M ItAR- ' »0< r 10 .5# : 
Ul^ X T f < A » X 2 I ) 


Of LOO I no 
bPLOO 1 VO 
CF L00200 
l»Pl 00210 
GF‘l 002 20 
(O L00230 
Gf I 00240 
D) 1.002S0 
Gf I.002A0 
Cil 1.00270 


121 FORMAT < //2*iX » ' C.ROUNH fATURN 

X > 

NUf ! s -;w 

FAl I PI OT ( A. #0. # - 3) 

Yl MfOX- KUF ' 0 
XM -XI MAX /U . 

x..( ’ 1 .0 

XM I N - O . 

1»XM IN-- 1 , 

i'O UOJ If'«l . JUO 

IF (XMiN.i I. MM on TO no2 

XM I N KM I N F »*XM I N 
UO 1 (UNIX NUI 
U02 XMAX-*Sl' 

;.X AX I S-»Ml XM IN 
NAP ^ S-IIXAX I *.♦ 1 

X AX I S ( I ) Xfi I N / XM 
1 AX I •*. < I > - O . 

IM- MOV 1 t -.••NAP I I, 

X AX j *»« II ) -xAx I s ^ II I > 1 1 . / x';l 
M^ K I ‘ * n I * > A X I S < If II 
(to ( ON I 1 NUt 

(MAX I / X'.( 

xA>i‘,(NAi >m: 

. AX 1 •. i NAl I ■ . ♦ I ' O . 

I»0 MOV 1 > ' ’ » NAl J i, 

X%X|*»<NAI T‘.MI )-Sl>^XS(. 
>.«Xi*.(N»W IMl. >.>AXl‘j<NA» I'lFZI 
no*. ( ONI lMt<( 

Nl { 01- 'tNAl I 

I At I t I Nl 1 v»iX M, f > AX X r> f Nl I (I T • 1 
i Al I •..YMI'I 4 ^ * . » tl • . • MM* - 1 4 

» I I I M . ✓ • -P. 

XI t Nil! I •. 1 


OKOUNK fATURN IN JF T (ENTERFD COOROINATF: SYS f E H ' / /Gf L002G0 

’ bPl 00290 

• 00-^00 

• - 3 ) Gl w A03 » 

r»r LOO i: o 
Gf I OOJJO 
GPl 00 T4<) 
(Vf L OOJM> 
1>I I 00 lAO 

. (,fnooT/o 

> on TO no2 c.iLoouio 

' GPI 00 JVO 

M I OO’i^O 

^ ni t 00 )io 

' (»l I C04.*0 

T.f 1.004 TO 

M (»f I 00440 

Gf I 004'. o 

I t, l.l U004 AO 

>\ll I > I 1 . /X';L Gf 1.004/0 

< n I I (.f i o04f;o 

Gl I 004VO 

(.1 ( OO'lOO 

'.;i» y *.i: I*f i oos I o 

o. (»l I f^OSJO 

1 t, l.l I JO*. 10 

-Sl>^ XS(, l.» I O0S40 

. > AX 1 ‘j < NAl )MZIMM|*VAX l.ll OO . .O 

* 1*1 I O0*-.A0 

(.1 ( 20*. /O 

'. » > AX xr. • Nl I 11 T • I . I • U I • . 1 > (.1 L 00 .MO 

. • tl • . • MM*-14 O. l./U=2.00 •0.»lin Gf I 

(.( l OOjsOt'* 
l.f [ Om • » <•> 
Gi 1 O0A.‘O 


U I • . 1 ) 

G/ G = 2 . 00 


Ii‘J,U 1 l' \-.t. C'O’.'tilUtL’ti 


ii-ir 


I iL SUkai,* 







|J»|* V 





DO ^OO i»l«HU 

URirr (Ar66) PU(1) 

66 FOKMMr <///4i,Xf 'PBARa' »riO.S) 

1FLA0« J 
XPUF«^0. 

IJc T*1 

lU-l 

RUaSl' 

ALFUG^t .5 

CALL UAI.LJ (Rb;»V>lN*0UH2*DU»3»DUn4f PEL PWi> 
PMAXO«PrLFW 

xr <F MAXG.L^ U< 1 > > XFLAC-O 
Ru*bp/c(ii; <rMUio> 

CALL UALL } ( K U » r Pi IH? f DUH3 f DUH4 # PEL PU > 

PrtAXF = prLI’U->^MN**7* ( VHNSl OS < PM 1 UP) )»»2 
IF <XrLAG»£0.0> GO TO SOL' 
k'^ALL ( 1 ) ^r>P 
r WALL < 2 > ^*>P* . 1 
DU 400 K«1 » lOO 

CAl L UALL J ( KWALL ( K ) • UMN» PUn2f DUM3f 0Uf14 » DLLF U > 
XU---SOK T < RWAi t < K ) *4? - spr«::» ) 
rHiu^AiAwroUf bp/ 

PM AX « PE LI U-UMN**2M VMN4C0S<PHIU> )*42 
FKK< K ) =-f Ut I > - PMAX 
IF ( K ♦ CO . X ' CO TCI 400 

IF <Apf:(CKR<K> ) .LT. I .t-b) GO TO 50X 
r>KL>ia(KUAt.L (K > 'RUALL<K-1 > > / < ERR ( K > > EKR < K - 1 ) > 
hUALL < K 4 1 ) •-KWAl L < K > 'REL AX»PRPE «CRR<K > 

400 rUNT ZNliL 
TUF 

SOI LU-RUAL L • K > 

\f U»“- t-G»x r ' KU**2 -SD» ) 
rur-AiANwMXi iir*sc») 

Yur-'cui 4MG1 *M»/LnS<PliF ) 

F Ul p'-f Uf e U4C . /r I 
SOJ I«l HiU f 111110/ XF TS 

ir <f HUio.LrKo. > GO u) joi 
f H I U O . 

r«j .'00 ; I* ; .mI'i ts 

0‘. . f H 1 U ^ 

C Al I WAl I J < F W # VMM » [*Lm J » Uim3 » PUM4 • DLLPW > 

YU < u 1 *M r. I • '^;p/cos «. r h i n > 


XU -M«* 1 AN <1 H I U ' 

I H I - A I AN.’ * <U * *'.fi - fU ) 

KU ( M. U > ✓ t Oil O HI > 

KUJ Ko 4 r I s* I . * rr)‘. ( PH I > > 

L TAM- KU 'KUI 

II ( f I AM (.1.1.) F»TAC.-0.0 

I r ( I 1 AM (, T 1 . ) (,0 , U I o T t 

4 ( I Al> I I I AU T AM . At I (. > 

1 O I 1 I M I N Ilf C I S4I L 1 At. 

r M.."* PFLI U OMN« ♦ 4 ( OMNttf ()<; I » H I U ) > *4 

(AM I MA U H < I M 1 N . f MA * . Al.4 UG > 

* l«AK - O U ( I ) r M 1 N ) / ( f MAX f M I N ) 

II 1 » H I N . G i . I U ( 1 > > l.U I U *.0-r 

. KA) J I U ( M , I'LLf b 
If ■ ( l-Af I . (, I I > GO UY .OV 

I I .. I k I f P tK J • • . J . ) ♦ • M . / AL 4 G ) 

* I i U 1 U 1 it T I KOI « I I A I* t ru. ( I H 

if ’ " • t . \ If I ) fvOf • i I A J ft ‘v 1 N < » MI » 

X I I m K 1 f 1 ' - X I t ,M I i I K t ) / • M 

f I U I 1 t 1 n I • r 4 I II 1 I M If I • ♦ xM 

I .» ? ' li 1 ♦ 1 

«>•' II ■ t !*AK . l. f ; (JK . I l«AK . I I O . * i.U 

‘ I A( I I 1 « WAK ♦ • . . , J * » i I /At I Ul. » 

vi ( u n M ( u « t u ♦ I I I 


0PL0O630 
GPL00640 
GPL00650 
GPL00660 
GPL00670 
GPLOOAOO 
GPL006V0 
GPL00700 
GPL00710 
GrU007Z0 
GPL00730 
GPL00740 
GFLO0/5O 
GPl 00760 
GPL0077O 
GPL00700 
GF L00790 
GPl 00800 
CPL00810 
GPL00020 
CPLOOG TO 
npl 0v'>84O 
OF LOOOSO 
GPL00060 
t»PLO0070 
GF I OOOGO 
CPLOOavO 
OPl 00900 
GF L00910 
GF L00920 
GF L00930 
GFLO094O 
GFL009S0 
CjF L00960 
GF L009/0 
GFL00980 
GPLOOVVO 
G‘^LOlOOO 
LOlOlO 
GFLOIOJO 
GPLO I 030 
0PL01040 
GPLOIOSO 
cr LO I 060 
GFL01070 
GFLOlOnO 
GPL01090 
GF LOl lOO 
OF LO 1 I I O 
GF LO I i .70 
Of LO I 1 50 
GF I 01 1 40 
GPL 01 1 GO 
GPLO I 1 60 
GF L 0 I 1 /»0 
GPl O ! 1 GO 
GF L O I 1 90 
C.F L O 1 TOO 
GPl 01 I'l ;> 

GF L o I rro 

GM O i .M J 
(.1 L O 1 ’4 > 
01 L O ) .'M> 

I L) 006 i.f I <• . .V.O 

(',» I 0 1 7/^0 

OF L o 1 ruo 


1 i^i’i I A-.;. C ontinuLil 


11-11 





fj ►J 


fa 


YFI.nTU< IU)»)iU 
XFLOTU( IU>=KFLOTU< IU>/XSC 
YPLoru< iu>*-Yr‘LOTU< un/xsc 
IU-»IUf 1 

PHIU«PHUUDPHXU 
CONT INUt 

HI ►tl = < P I - f HlO> / XF TS 
I’MIeF HIO 

DO AOO IK*'l*NUFTS 

rhfiifi 1 >/un.f s 

IF O l»T . i . ) CO TO AOO 

t T*4 C 1 • - I L<AKJ«« . 25 >«a ( 1 . /AuPC / 
Ror-‘ROftr3»/ r 1 . ♦coscFHi ) ) 

XPLOTJ< IJCl ) »KOF OlCTAjacaScrMI ) 

Vf LOT J< IJC r > -ROPaE TA J»S I Ml PHI ) 

XF I or j< i jr r > “XF LOT j< ijE T ) / X5C 
YF I OTj( uei )-»Yf LOT J< IJE 1 >/XSC 
iJt T ^ ur. T ♦ 1 
r MI Ml f ItF'H I 
PONT INUt 

IF < in Ab.LO.O) GU TO 513 

IF <ru( 1 ) .EO.O. > G«1 TO 513 

IF < r OF . I T .1 MlUO ) CO in 510 

GO TO 51 t 

xf L onii un 

Yf l OUi< lin-YFUF 

XPL OUl < 1 U > *• XPL (»t LI < lU ) / XSC 

Yr>lO?U<IU)«YMOru<IU)/XSC 

no TO uoo 

r»F M TO - < MIF - f HZUO ) / XPTS 
I HIU--I HUH) 

no ^00 ij- 1 pNut Tb 
Yu-coi *‘j 10 1 ai,n/('os<PHiU) 

XU-SnaTAN<» MIU) 

Ku~-‘*n 'c'<^o<n<u«> 

( At 1 UAL I J i RU* VMN * 0UM2 » 00^3 r IHJM4 t L'EL f U > 
f MAx«.ru,i pu-Of-Naa.M <OMN*cas<PHiu) >a*2 

f I.AF - lU \ 1 ) ^ I HAX 

L" TAO ( 1 . -p HAL a * . .5 ) ** < I . /ALPUG ) 

XF t OTiu lu) *5I»' roat I AtJ 
Y n u I u< I u » - XU 
xF \ o 1 u( I II ) - xf 1 n nj < no / xsr 
YF LOT U < lU ) •» VF LOTH ( I o ) / XSC 
lU- lUF 1 

F‘HIO-FmUUIiF HIU 
( TINT I NUL 
I U - I U I 
1 .Jl 1 -- IJL J - I 


Of LOl .'90 
(»( I 01 TOO 
CF LOl BIO 
Cfl LOl 320 
GPLOl 330 
CPL01340 
GPLOl 350 
GPLOl 3AO 
GPLOl B70 
GPLOl 3HO 
GPLOl 390 
GPLO 1400 
GPLOl 4 I O 
GPLOl 420 
r.PLOX 430 

Of LOl 4 40 
DF'LO I 450 
nr LOl 4 AO 
CF LO 1 4 /O 
GF LOl 41)0 
GPLOl 490 
CPI. 011.00 
GPLOIMO 
Gf I 01520 
GPl 015 TO 
CF LO1540 
GPLOl 550 
GPLOl •'.AO 
bl'L 015/0 
CF LOl'iHO 
GF I 01*. VO 
GF L O 1 AGO 
GFLOt AlO 
GPLO I A20 
GPl OS A30 
GF I (; 1 A40 
GF I O I A50 
GF I Oi AAO 
GfLOlA O 
GPl OldllO 
GF I O I AVO 
C>F I O I ’OO 
GF L O 1 / I O 
GF LOl /20 
IW LOl / Bo 
GF I O I / 40 
GF I O 1 /50 
Cri.O I /AO 


UK I rr 

< A . 1 ) > 1 

It. T nij 



GF 


Oi //o 

t 1 F OF 

MA 1 

< 1 5X . 

JfT IMFINOfMENT REGION* 

IJET-*'*I4#10X# 'Uf WA*»H 

nrrLijF 


O 1 /DO 

U t T 

I ON 

M G 1 lIK . 

HI '.14//) 



GF 


O 1 /VO 

I f 

( 1 FF 

7.GI.O.OK.UI G1.0' UF^lTP < At 552) 



GF 


O 1 FIOO 

5'* F tlF^f-A? 

l /• / 3 ♦ X » 

X150.) .AX# YiGU.i . nx . ' X isnin #ox# fiGOU /> 


GF 


O 1 t> 1 o 

I n I N ' 1 u 




GF 


o I H .*0 

1 l 

IF 1 




(»F 


0 1 D BO 

1 1 

V U) 

(.1.1 IF T 

> IMAX-UI 



GF 


O 1 D40 

X F 

( U 

. 1. 1 . 1 It 1 

) I H T N - I .IE T 



(>l 


0 1 D50 

no 

•lO B 

1 M - 1 . J “ 

\X 



(.F 


C> 1 t • /. O 

1 F 

< M 

.1 t . IMIN 

» UM Tl ( A# 1 2 ; XPL OT 1 < IM ) . 

YPl OT J< IM > # Xf LOl U ( 

in > # YPLon'.F 


O I D 'O 

l 11 • t 

M 





GF 


O 1 DDO 

1 2 F UFx 

MA 1 

( BO' . 4 * 

1 O . 5 . 2X ) ) 



GF 


O 1 D*. O 

1 F 

( M 

. l> I . 1 M I R 

. »V.n . I M 1 N . F O • 1 U ) UF^ I 1 t <A 

• 15 ) XF'l OI M I « ) # Y 1 

1 Ot J< 

in) GF 


0 1 /Oo 

1 5 F OKnA I 

. BO X . . ( F I 0 . , 2X ) ) 



(G 


O 1 •’ 1 V) 

I F 

4 1 W 

(.1 . 1 M 1 R 

. ANI» . 1 f-. 1 N.KJ . 1 JE 1 ' UKIIE 

( A • t A > Xf t oni< 1 n > . 

rl 1 O’ 

U( im.i 


O 1 V 2 o 


(.» I vM V 
L.F I O 1 V >0 


} r;uio A-.^. 1. L’r*n'u<‘(i 


li-12 





e03 CCiNTtNIiC: 

C CALI I INF ( XF I 01 Uf YF LOT U » 111* 1 » 1 f i * I » « 1 > 

C CALL LTNL < XT L Q 7 J # Y< L O T J • ZJC T • 1 » 1 • * • t f . 1 ) 

SOO LON7INUL 

Vt NIO- . Off I/XPI iJ 
PMIU^O. 

un voo I i •« I #NlTf IS 
YU^CIII^SIGT HHJ) 

XO=-5»i»® 1 AN .7 ill U> 
ir < aU . liT . Xl.liAX > GO TO 901 
XI LO UK I J) *bl« YU 
Y» I LUIK 1 D-XU 
XllOIlKI I > - XI I IJ I U< I J ) /\GC 
YI‘LO UK I J > V. YT'LO TU< IJ > / XGC 
I tU (i -F M I Ui t'F HIU 
VOO rONTINUL 
VO 4 I IJ- I J I 

UK 1 U < (lO / > 

U07 f OKNAT <,//lSX* Ul UASH ULrCCCIlON 70NL LINE# fl«AR-«0 OUlSIl'f 
ici ION kh.ion // > 
ut. 1 1 1 ( A • lion > 

OOU MJKMAT 1 ibX . ' XUP ' * HX • ' I'UP ' // ) 
lUl U I IN- I . lU 

URITC X> L OliK IN } * > f'L O UM 2n > 

n * ^ LON } I NLU 

til I I OKMA t i SOX . .U r * O. 5 I ?X ) I 

( Ai I LINT • <1 LU t Ur Yf (. OT U« ll/« 1 « 1 » t « 1 • . I ) 

(' (AM A(»KAW 

L I.Al L 1107 <-V9.r VV.r~T> • 

r Ki AO ( br lOl ) F AU'jt 
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